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ABSTRACT 


Floodplain  encroachment  has  many  ramifications.  The 
present  study  investigates  the  possibilities  of  scour  and 
excessive  backwater  resulting  from  the  construction  of  the 
Columbia  Generating  Station  on  the  Wisconsin  River  flood- 
plain  near  Portage,  Wisconsin.  Neither  scour  nor  excessive 
backwater  appear  to  be  problems  in  the  study  area,  not  even 
for  the  100-year  (regional)  flood.  A possible  change  in 
climate  since  1950  does  not  alter  the  expected  results  of 
floodplain  encroachment.  However,  the  effects  of  the  pro- 
posed climatic  change,  manifested  by  a shift  from  predomi- 
nantly zonal  to  predominantly  meridional  circulation  pat- 
terns, cannot  be  fully  determined.  Limited  data,  few 
choices  of  analytical  methods,  and  the  possible  masking  of 
any  effects  by  reservoir  development  all  contribute  to  the 
inconclusiveness  of  the  analysis  of  climate  change  and  its 
possible  effects  on  the  study  problems. 
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INTRODUCTION 

This  thesis  analyzes  problems  of  increased  water  surface 
elevation  and  increased  flow  velocity.  The  problems,  arising 
during  flood  flow,  result  from  man's  encroachment  on  a flood- 
plain.  Increased  water  surface  elevation,  a result  of  the 
backwater  effect,  could  possibly  cause  overtopping  of  flood- 
plain  structures.  High  stream  velocities  increase  the 
river's  potential  to  scour  the  floodplain  surface  and  fea- 
tures. Study  of  these  hydraulic  problems  initially  takes  ac- 
count of  current  climatic  and  hydrologic  data.  However, 
increasing  evidence  shows  that  Wisconsin,  like  most  of  the 
world  for  that  matter,  is  being  affected  by  a change  in  cli- 
matic regime  (Lamb,  1966;  Kalnicky,  1974;  Knox  et  al . , 1975). 
In  an  attempt  to  evaluate  this  possibility,  climatic  and  hy- 
drologic records  for  the  periods  1873-1910  and  1925-1975  are 
examined  for  evidence  of  climatic  change.  The  two  geomorphic 
problems  are  re-examined,  as  hydrologic  characteristics  rep- 
resentative of  the  actual  climate  regime  are  considered  in 
the  reappraisal.  The  results  of  the  present  research  demon- 
strate two  major  themes:  the  effects  of  man  on  a stream's 
flood  flow  regimen,  and  the  effects  of  climatic  change  on  the 
study  problems. 

ORIGIN  OF  THE  PROBLEM 

Floodplains  can  be  troublesome  areas  for  people,  espe- 
cially floodplain  managers,  planners,  and  inhabitants. 
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Conversely,  floodplains  have  historically  been  optimal  places 
to  live,  whether  because  of  proximity  to  water,  of  transpor- 
tation routes,  or  of  the  location  of  factories.  Floods  can 
more  than  counteract  the  benefits,  however.  Kates  (1962) 
investigated  the  reasons  why  people  lived  on  floodplains  and 
their  perceptions  of  safety  versus  danger  and  benefits  versus 
detriments.  Apparently,  people  are  willing  to  make  trade- 
offs. They  either  accept  or  ignore  the  dangerous  possibili- 
ties in  orde..  to  receive  the  benefits  of  living  on  a flood- 
plain.  Thus,  floodplain  development  and  occupation  continue 
to  increase.  Problems  are  compounded  when  one  considers  that 
floodplain  occupation  often  adversely  affects  the  hydrologic 
regimen  of  rivers.  Two  results,  according  to  Wolman  (1967), 
are  increased  runoff  and  sediment  yield  (especially  during 
the  construction  phase) . In  an  attempt  to  counter  indiscrimi- 
nate floodplain  development,  state  floodplain  regulations  and 
a National  Flood  Insurance  Program  have  been  established  in 
the  United  States.  Unfortunately,  floodplain  zoning  ordi- 
nances are  often  vague  and  difficult  to  enforce.  Enforcement 
difficulty  is  manifested  in  the  situation  which  occurred  at 
the  study  site  of  the  present  research. 

The  Columbia  Generating  Station,  a coal-fired  power 
plant,  began  operation  in  March,  1975.  This  527-megawatt 
station  is  located  on  the  Wisconsin  River  floodplain,  approxi- 
mately 6.4  km  south  of  Portage,  Wisconsin,  in  Columbia  County 
(See  Figure  1) . There  was  much  opposition  to  the  construction 
of  this  plant  when  it  was  first  proposed  to  the  public  in  1971. 
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Environmentalists  and  other  concerned  parties  took  advantage 
of  the  then  recently  approved  federal  law,  the  National  Envi- 
ronmental  Policy  Act  (NEPA) , in  order  to  express  their  mis- 
givings. Most  were  concerned  with  the  possible  atmospheric 
and  biologic  pollution  which  could  result  from  the  sulfur  di- 
oxide (SO^)  content  in  the  coal  to  be  used  in  energy  produc- 
tion. 

However,  others  were  concerned  because  the  plant  would 
by  wholly  situated  on  part  of  the  Wisconsin  River  floodplain. 
Thomas  Lee,  a hydrologic  engineer  with  the  Wisconsin  Depart- 
ment  of  Natural  Resources  (DNR) , voiced  his  opposition  to  the 
construction  of  the  plant  at  a 1971  public  hearing.  In  his 
opinion,  the  plant  was  in  violation  of  a 1965  DNR  floodplain 
zoning  ordinance,  in  that  it  would  create  a greater  than  al- 
lowed backwater  effect  (Schmied,  1973).  The  Environmental 
Impact  Reports  and  Statements  associated  with  the  plant  con- 
struction also  alluded  to  backwater  as  a potential  problem, 
in  addition  to  other  hydrologic  phenomena,  such  as  increased 
flow  velocities,  flow  pattern  changes,  scour  and  fill  pattern 
changes,  and  possible  scour  of  the  plant's  cooling  lake  dike 
(EIS,  1974). 

The  various  objections  and  doubts  went  unheeded  as  the 
Columbia  County  Board  of  Adjustments  granted  a variance  to 
the  floodplain  ordinance  in  favor  of"' the  power  plant.  Possi- 
bly, the  tax  island  which  could  be  established  by  the  utility 
tax  revenues  meant  more  to  the  officials  than  did  their  en- 
vironment.  The  impact  of  the  tax  benefit  was  studied  by 
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Michael  Shaver,  a former  student  in  the  University  of  Wiscon- 
sin's Geography  Department  (Shaver,  1976).  Whatever  may  be 
the  reasons  behind  the  decisions  involved,  the  final  result 
was  the  appearance  of  the  Columbia  I Generating  Station  on 
the  floodplain  of  the  Wisconsin  River. 

In  the  final  analysis,  it  seems  that  the  progress  of 
man  must  prevail  over  the  preservation  of  the  natural  envi- 
ronment. All  too  often,  however,  alternative  solutions  to  a 
problem  are  either  not  considered  or  are  ignored.  Did  the 
Columbia  Generating  Station  have  to  be  built  on  the  flood- 
plain?  Could  it  have  been  situated  farther  away  from  the 
river,  off  the  floodplain,  and  still  accomplish  its  purpose? 
These  questions  are  now  moot  in  the  case  of  the  Columbia 
power  plant,  but  they  can  and  should  be  carefully  considered 
in  future  developments. 

The  present  research,  in  addition  to  the  backwater 
problem  raised  by  Mr.  Lee,  examines  one  of  the  major  problems 
discussed  in  the  Environmental  Impact  Statement  (EIS) : the 
effect  of  increased  flow  velocities.  These  hydraulic  prob- 
lems are  considered  not  only  in  relation  to  the  study  site, 
but  also  with  a view  to  developing  a method  of  analysis  ap- 
plicable to  other  sites  of  the  same  general  kind. 

Climate  adds  another  dimension  to  the  discussion.  Knox 
and  others  (1975)  demonstrated  for  the  Upper  Mississippi  River 
Valley  region  that  there  is  a significant  persistence  in 
climatic  regimes.  Their  study  refined  the  work  done  by  Lamb 
(1966)  on  a hemispheric  scale  and  by  Kalnicky  (1974)  on  a 
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continental  scale.  All  proposed  the  need  to  consider  climate 
as  being  essentially  "nonrandom".  Lamb  (1966)  postulated 
that  shifts  in  climatic  regime  are  due  to  a change  in  the 
tropospheric  wind  patterns.  A shift  is  said  to  have  occurred 
around  1960,  as  the  long  wave  patterns  shifted  from  a pre- 
dominantly zonal  (westerly)  flow  to  a more  meridional  (north- 
erly and  southerly)  flow.  The  meridional  regime  apparently 
resembles  the  cooler  and  moister  pattern  of  general  circula- 
tion that  existed  prior  to  1895.  The  results  of  dominant 
meridional  circulation  are  persistent  periods  of  high  year- 
to-year  variability  in  temperature,  precipitation,  and  hy- 
drologic response,  such  as  larger  floods  and  sediment  yields 
(Knox  et  al . , 1975).  For  the  Upper  Mississippi  Valley  region, 
droughts  and  extreme  floods  are  alike  characteristics  of  me- 
ridional regimes,  whereas  the  frequent  recurrence  of  droughty 
weather  dominates  in  predominantly  zonal  regimes. 

Significantly,  this  recent  shift  of  climatic  regime  pos- 
sibly establishes  the  facts  that  climate  change  is  episodic 
and  that  weather-related  events  are  nonrandom.  Since  most 
historical  rainfall  and  streamflow  records  date  back  only 
100  years  at  best,  it  is  apparent  that  most  of  the  data  are 
biased  toward  the  zonal  flow  regime  that  existed  from  about 
1895  to  1960.  Because  of  this  zonal  domination,  the  30-year 
climatic  normal,  1930-1960,  may  be  the  most  abnormal  "normal" 
that  could  have  been  used  for  the  aggregation  of  climatic  and 
hydrologic  data.  A question  must  arise:  how  representative 
are  the  existing  hydrologic  data  and  determinations  of  flood 
recurrence  intervals  used  by  floodplain  managers,  floodplain 
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construction  engineers,  and  designers  of  flood  control  struc- 
tures? The  present  research  attempts  to  develop  the  topic, 
posed  by  the  above  question,  in  the  examination  of  the  two 
geomorphic  problems  discussed  previously. 

Both  problems,  backwater  and  increased  flow  velocities, 
may  occur  because  of  the  power  plant’s  occupation  of  space 
formerly  used  for  the  transport  and  storage  of  floodwaters. 
The  plant's  presence  on  the  floodplain  creates  a constriction 
which,  during  high  water  periods,  causes  increased  flow  ve- 
locities near  the  constriction  and  increased  water  surface 
elevations  upstream  and  at  the  constriction.  The  increased 
velocities  could  possibly  cause  scour,  which  could  damage  the 
floodplain  surface  and  structures,  lead  to  increased  sediment 
load,  and  create  ecological  hazards,  such  as  sediment  deposi- 
tion on  fish-spawning  grounds  downstream.  Increased  water 
surface  elevation,  the  result  of  too  much  backwater,  could 
cause  overtopping  of  upstream  structures.  Though  other  sig- 
nificant potential  problems  were  raised  by  the  Corps  of  En- 
gineers in  their  EIS , the  present  research  is  limited  to  the 
problems  aforementioned.  This  researcher  believes  that  the 
study  problems  clearly  demonstrate  bow  man  can  alter  the  hy- 
drologic regime  of  a river  during  floods.  The  two  problems 
under  investigation  apply  to  any  construction  on  a river 
floodplain.  Since  many  rivers  attain  bankfull  stage  about 
once  a year  on  the  average,  and  inundate  their  floodplains 
at  quite  modest  recurrence  intervals,  the  phenomena  created 
by  the  power  plant's  presence  can  recur  at  many  other  sites. 
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The  present  research  fills  a void  in  the  work  of  the 
Columbia  Site  Impact  Study  Group.  This  organization,  under 
the  leadership  of  Dr.  Daniel  Willard,  was  formed  to  study  the 
environmental  impacts  resulting  from  the  construction  and 
operation  of  the  power  plant.  Since  his  group  began  the 
study  in  1971,  Dr.  Willard  has  expressed  a desire  for  an  ex- 
amination of  the  changes  in  surface  water  flow.  The  present 
study  attempts  that  examination. 


REVIEW  OF  THE  LITERATURE 

Three  major  topics  are  considered  in  the  present  re- 
search: scour,  backwater  and  climate  change.  For  purposes  of 
this  discussion,  scour  means  entrainment  and  movement  of 
particles  as  a result  of  the  flow  of  water  across  them. 
Backwater  is  the  increase  in  the  water  surface  elevation  of 
a river  due  to  the  effects  of  confluences  and  constrictions 
on  the  flow  of  floodwaters.  Climate  change  is  a shift  in  the 
dominant  upper  air  circulation  patterns  with  resultant  changes 
in  weather-related  events. 

The  Columbia  Generating  Station  forms  a constriction  on 
the  Wisconsin  River  floodplain  during  periods  of  high  water. 
Flow  velocity  (V)  must  increase  at  the  constriction,  because 
the  presence  of  the  power  plant  reduces  the  cross-sectional 
flow  area  (A)  for  a given  amount  of  water  discharge  (Q) 

(the  statement  follows  from  the  fact  that  Q=AV) . Increased 
flow  velocities  provide  increased  erosive  ability  for  the 
flowing  water.  Therefore,  the  floodplain  surface  and  fea- 
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tures  near  the  constriction  should  be  more  liable  to  scour 
than  are  other  areas  nearby,  assuming  that  soil  particle  char- 
acteristics are  sensibly  uniform  throughout  the  area  in  ques- 
tion. Scour  depends  upon  flow  hydraulics,  as  well  as  on  the 
particle  size,  shape,  density,  and  cohesiveness  (Bryan,  1976). 

Flow  hydraulics  have  been  expressed  by  Leliavsky  (1955) , who 
indicates  three  controls  on  scour:  critical  flow  velocity, 
critical  tractive  force,  and  the  buoyant  or  lift  force.  Both 
Leliavsky  and  Graf  (1971)  agree  that  the  field  investigation 
of  the  lift  force  is  impracticable  at  the  present  time. 

Critical  flow  velocity  seems  to  be  the  most  useful  criterion 
to  consider,  since  it  is  fairly  easy  to  employ.  However,  un- 
certainties exist.  There  is  a question  as  to  what  flow  ve- 
locity should  be  used.  Some  researchers  have  used  the  bottom 
flow  velocity,  while  others  have  employed  the  mean  velocity 
as  their  threshold  variable  (Novak,  1973).  Graf  (1971)  also 
cautions  against  indiscriminate  use  of  critical  velocity  charts, 
such  as  that  of  Hjulstrom  (1935).  If  study  conditions  re- 
semble the  empirical  conditions  used  by  researchers  into 
critical  velocity,  then  the  method  may  be  employed.  Since 
the  present  research  entails  similar  conditions,  such  as  non- 
cohesive  soils  and  comparable  particle  sizes,  the  critical 
velocity  threshold  has  been  selected  for  investigation  of  the 
scour  potential  of  the  floodplain  surface  and  features. 

Floodplain  scour  has  been  noted  in  many  studies  (Jahns, 

1947;  Schmudde,  1963;  Nelson,  1965).  In  all  cases,  however, 
the  extent  of  scour  was  small,  being  limited  to  local  occur- 
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rences.  A low-frequency,  high-magnitude  flood  event  (50-100 


year  flood)  did  not  appreciably  scour  the  floodplain  (Gupta 
and  Fox,  1974).  Dury  (1973)  found  the  same  lack  of  scour  for 
an  even  more  extreme  event- -a  flood  with  a recurrence  inter- 
val of  between  500  and  1000  years.  Since  the  possibility  of 
floodplain  scour  appears  to  be  limited  to  small  and  local  oc- 
currences, map  and  visual  reconnaissances  near  the  power  plant 
establish  the  area  near  cross-section  8 to  be  the  most  sus- 
ceptible to  scour  (See  Figure  2) . Here  the  normal  flow  of 
flood  waters  is  greatly  constricted  by  the  west  dike  of  the 
power  plant's  cooling  lake  and  a small  knoll  west  of  the  dike. 
The  constriction  near  cross-section  8 creates  a narrow  flood 
channel.  Scour  should  occur  here,  if  it  is  to  occur  at  all. 

For  that  reason,  the  investigation  of  the  possibility  of 
scour  of  the  floodplain  surface  and  features  is  limited  to 
the  area  that  includes  cross-section  8. 

Although  the  backwater  effect  is  a well-known  hydraulic 
phenomenon,  it  is  treated  in  only  a cursory  manner  in  the 
engineering  textbooks  of  the  West  (Linsley,  Kohler,  and  Paul- 
hus , 1958;  Ward,  1967).  Considering  the  variety  of  studies 
on  record  (Rozenberg,  1972;  Walik,  1971;  Lodina  and  Chalov, 
1971) , Soviet  hydrologists  seem  to  regard  backwater  problems 
to  be  important  areas  of  research.  In  engineering  practice 
at  the  state  and  federal  levels  in  the  United  States,  com- 
puter programs  have  been  written  to  calculate  the  backwater 
changes  resulting  from  encroachments  on  the  river's  flow  area. 
For  example,  the  state  of  Iowa  has  its  own  program;  Wisconsin's 
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DNR  uses  the  Corps  of  Engineer's  model;  while  the  U.S.  Geo- 
logical Survey  has  yet  another  program.  All  programs  do  basi- 
cally the  same  thing.  Bernoulli's  Theorem  and  the  Manning 
Equation  are  used  to  compute  flow  changes  between  successive 
cross-sections  in  a river  reach. 

Many  state  floodplain  ordinances  limit  the  amount  of 
backwater  that  may  be  produced  by  a new  floodplain  develop- 
ment. Wisconsin  places  a 15.24  cm  limit  on  the  increase  in 
water  surface  elevation  near  urban  areas  for  the  regional 
(100  year)  flood  (See  Figure  3).  Mr.  Lee  of  DNR  calculated 
a 15.54  cm  increase  upstream  near  the  city  of  Portage,  be- 
cause of  the  encroachment  of  the  Columbia  power  plant.  How- 
ever, the  Corps  of  Engineers  and  the  power  plant  construction 
engineers,  Sargent  and  Lundy,  arrived  at  lesser  figures. 

Mr.  Lee's  testimony  was  not  accepted  by  the  Columbia  County 
Board  of  Adjustments. 

Dr.  Willard  of  the  Columbia  Site  Impact  Study  Group  has 
been  concerned  with  the  disagreement  about  the  backwater  ef- 
fect since  the  1971  public  hearings.  The  present  research  re- 
examines the  question.  The  Corps  of  Engineer's  backwater  sur- 
face profile  computer  model,  HEC-2,  is  the  research  tool  used 
for  the  analysis.  This  program  is  employed  by  the  Wisconsin 
DNR  today,  replacing  the  state  of  Iowa's  program  used  by 
Mr.  Lee  in  his  analysis. 

Climatic  change  is  a controversial  subject.  Abundant 
theories  attempt  to  relate  long-term,  short-term,  periodic, 
and  episodic  changes  in  climatic  regimes  and  weather-related 
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events  to  some  variable  or  group  of  variables.  The  present 
research  is  concerned  with  a recent  change  in  climatic  regime. 
Since  about  1960,  many  researchers  have  noted  a general 
cooling  trend  in  the  northern  hemisphere.  Surface  and  ocean 
temperature  anomalies  and  an  increase  in  the  southern  extent 
of  the  polar  ice  cap  support  the  belief  in  climatic  change 
(Lamb,  1972).  Earlier,  Lamb  (1966)  proposed  a shift  in  the 
predominant  tropospheric  long  wave  patterns  as  the  reason  for 
this  change  in  climate.  These  upper  air  circulation  patterns 
are  the  driving  forces  behind  the  climate  at  the  surface  of 
the  earth.  Lamb  propounds  that  meridional  components  in  the 
long  wave  patterns  are  now  the  predominant  situation.  These 
northerly  and  southerly  components  account  for  the  high  year- 
to-year  variability  in  temperature  and  precipitation  which 
has  characterized  the  past  two  decades.  The  prior  climatic 
regime  was  predominantly  zonal.  Thus,  westerly  winds  domi- 
nated the  period  prior  to  about  1960.  In  North  America, 
westerly  winds  contribute  to  the  orographic  precipitation 
along  the  West  Coast  of  the  United  States  and  a generally 
warm,  dry  climate  throughout  the  rest  of  the  country,  espe- 
cially in  the  Great  Plains. 

Lamb's  hypothesis,  that  a shift  to  a predominantly  me- 
ridional regime  accounts  for  the  generally  cooler  and  wetter 
climatic  conditions  since  about  1960,  has  proponents  and  op- 
ponents alike.  Some  opponents,  such  as  Willett  (1975),  do 
not  even  accept  the  dominance  of  meridional  components  in  the 
tropospheric  wave  patterns,  let  alone  Lamb's  hypothesis. 
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Recent  data  cloud  the  issue  somewhat.  Brinkmann  (1976)  found 


that  temperatures  in  North  America  have  increased  0.2°C  in 
the  past  decade.  Mason's  speech  to  England's  Royal  Meteoro- 
logical Society  establishes  the  facts  that  the  polar  ice  is 
now  retreating  northward  and  that  the  colder  temperature 
anomalies  are  decreasing  in  frequency  in  the  high  latitudes 
of  the  northern  hemisphere  (Mason,  1976).  However,  he  does 
concede  that  the  cooler  and  wetter  trend,  although  of  lesser 
intensity,  continues  in  the  mid-latitudes.  Since  the  upper 
air  patterns  control  surface  weather  in  the  United  States, 
Lamb's  hypothesis  is  a reasonable  statement  for  the  tempera- 
ture and  precipitation  anomalies  recorded  in  the  United  States 
and  investigated  by  Kalnicky  (1974),  Namias  (1970),  and  Knox 
and  others  (1975) . A noticeable  increase  in  meridional  long 
wave  patterns  has  occurred  in  the  Upper  Mississippi  River 
Valley  area  since  about  1950  (Knox  et  al. , 1975).  The  same 
situation  is  assumed  to  hold  true  for  the  Upper  Wisconsin 
River  Valley,  the  study  area  for  the  present  research,  be- 
cause of  its  proximity  to  the  Upper  Mississippi.  With  this 
assumption  as  background,  the  present  study  investigates  one 
of  the  climatic  variables  which  characterizes  a meridional 
regime--increased  precipitation  characterized  by  high  year- 
to-year  variability.  Kalnicky  (1974)  and  Namias  (1970)  fo- 
cused their  studies  on  temperature  anomalies,  whereas  Knox 

and  others  (1975)  examined  precipitation  characteristics. 

Lamb  (1972)  states  that  precipitation,  better  than  tempera- 
ture, reflects  the  characteristics  of  the  tropospheric  wave 
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patterns.  Thus,  only  precipitation  data  for  the  Upper  Wis- 
consin River  Valley  are  examined  in  the  present  investigation. 
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CHAPTER  II 

THE  STUDY  AREA 

Although  the  Columbia  Generating  Station  is  the  primary 
area  of  concern  in  the  present  study,  the  entire  Wisconsin 
River  drainage  basin  above  the  power  plant  must  be  considered. 
'A  downstream  site  investigation  should  not  be  undertaken  in 
isolation  from  a consideration  of  upstream  conditions'  (Cooke 
and  Doornkamp,  1974,  p.  88). 

THE  UPPER  WISCONSIN  RIVER  VALLEY 

The  Wisconsin  River  rises  in  the  Lac  Vieux  Desert  along 
the  Wisconsin-Michigan  border  (See  Figure  4).  Thence  the 
river  flows  southward  for  about  482  km  to  Portage,  then  turns 
westward  and  flows  approximately  190  km  to  the  Mississippi 
River  at  Prairie  du  Chien,  Wisconsin.  Elevation  decreases 
from  502.9  m at  the  river's  head  to  184.1  m at  its  mouth, 
with  greatest  rate  of  slope  decrease  in  the  upper  third  of 
the  river's  length.  The  river  flows  through  glaciated  and 
non-glaciated  parts  of  Wisconsin;  however,  the  reach  from  Lac 
Vieux  Desert  to  the  power  plant  is  situated  wholly  within 
glaciated  terrain. 

Above  Nekoosa,  resistant,  crystalline  rocks  of  Precam- 
brian  age  form  numerous  rapids.  Softer  sandstone  of  the  Upper 
Cambrian  Group,  covered  with  sand  and  gravel,  account  for 
shifting  sandbars  throughout  the  reach  from  Nekoosa  to  the 
Columbia  plant  (Herron,  1917;  Mead,  1911). 
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Level  and  undulating  topography,  except  for  a few  no- 
ticeable ridges  such  as  the  Baraboo  Hills,  prevails  in  the 
drainage  basin.  River  bluffs  are  few,  but  those  at  the  Wis- 
consin Dells  and  the  narrows  at  Dekorra  are  quite  spectacular. 

Lawrence  Martin  (1932)  divided  the  drainage  system  of 
the  Upper  Wisconsin  River  into  two  parts,  with  the  city  of 
Merrill  as  the  point  of  division.  Many  lakes  and  swamps  and 
a drainage  network  described  as  "crooked,  systemless  stream 
courses"  characterize  the  area  north  of  Merrill  (Martin,  1932) . 
A systematic,  dendritic  drainage  pattern  with  no  lakes  des- 
cribes most  of  the  basin  south  of  Merrill.  However,  from 
Wisconsin  Rapids  to  Portage,  an  "aimless  pattern"  with  large, 
undrained  interstream  areas,  few  tributaries,  and  a shallow 
valley,  replaces  the  dendritic  pattern.  Martin  (1932)  hy- 
pothesized that  the  absence  of  tributaries  and  a well- 
defined  valley  may  be  due  to  the  short  time  during  which  the 
river  has  occupied  this  part  of  its  course.  The  drainage 
near  Portage  is  thought  to  be  late  glacial  in  age,  as  a re- 
sult of  a change  in  the  Wisconsin  River's  course  from  the 
Devil's  Lake  gorge  to  its  present  position. 

The  river  water  is  reddish-brown  in  color  because  of  in- 
dustrial by-products  and  inorganic  pollutants  (EIS,  1974). 

The  sediment  load  of  the  Wisconsin  River  is  thought  to  be 
composed  mostly  of  bedload.  Though  no  bedload  data  are  avail- 
able for  the  river,  qualitative  estimates  suggest  that  sig- 
nificant quantities  of  sand  are  transported--enough  to  cause 
troublesome  recreational  and  developmental  problems  in  many 
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areas  (UMCBS,  1970).  Suspended  sediment  data  for  some  of  the 
Wisconsin's  tributaries  show  that  the  suspended  load  is  very 
low,  much  below  state  and  national  averages.  From  Lac  Vieux 
Desert  to  Wausau,  the  suspended  sediment  yield  is  only  3.5 
tonnes/km^  per  year,  while  the  load  increases  to  10.5  tonnes/ 

O 

km  per  year  in  the  stretch  from  Wausau  to  Portage.  These 
suspended  load  figures  are  well  below  the  state  average  of 
28  tonnes/km^  per  year  (Hindall,  1976).  Apparently,  sus- 
pended sediment  load  is  minimal  in  glaciated  parts  of  Wiscon- 
sin, but  often  significant  in  the  non-glaciated  Driftless 
Area  of  the  southwestern  part  of  the  state.  No  analyses  of 
particle  size  distribution  for  suspended  load  of  the  Wiscon- 
sin River  are  available,  but  fine  clay  sizes  probably  dominate. 
Small  percentages  of  fine  sand  and  silt  should  also  be  pres- 
ent (UMCBS,  1970). 

The  average  annual  runoff  decreases  from  33.0  cm  in  the 
north  to  20.3  cm  near  Prairie  du  Chien  (UMCBS,  1970).  Most 
flood  flows  on  the  Wisconsin  occur  in  the  spring  as  a result 
of  snowmelt  runoff.  However,  the  disastrous  floods  on  record 
have  occurred  in  the  late  summer  and  early  fall  because  of 
large  amounts  of  rainfall. 

The  flow  of  the  Wisconsin  River  is  heavily  regulated. 

Above  Portage,  there  are  21  reservoirs  (493,200,000  m^  stor- 
age capacity)  and  3 power  dams  (328,594,500  storage  ca- 
pacity) . These  24  reservoirs  and  flowages  stabilize  river 
discharge  for  purposes  of  power  utilization  and  recreation. 
Additionally,  many  small  hydroelectric  dams  and  3 large  flow- 
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ages  (at  Petenwell,  Castle  Rock,  and  Wisconsin  Dells)  exist 
north  of  Portage  (UMCBS,  1970).  The  reservoirs  and  flowages 
provide  a significant  amount  of  flood  flow  regulation  in  the 
spring.  Of  great  importance,  however,  is  the  fact  that  these 
dams  and  reservoirs  have  virtually  no  effect  on  flood  flows 
at  Portage  during  the  summer  and  fall,  because  at  that  time 
the  reservoirs  and  flowages  are  full  (Sherrar,  1976a). 

THE  COLUMBIA  GENERATING  STATION 
The  power  plant  is  6.4  km  south  of  Portage,  in  the  ex- 
treme southeastern  part  of  the  Central  Sand  Plain  (Martin, 
1932).  A unique  drainage  situation  exists  at  Portage,  where 
a distance  of  only  2.4  km  separates  the  Wisconsin  and  Fox 
Rivers  (See  Figure  1).  During  normal  flow  periods,  the  Wis- 
consin River  is  1.83  m higher  than  the  Fox;  the  height  dif- 
ference may  increase  up  to  6.10  m during  extreme  flood  events 
(C  of  E,  1971).  Historically,  interbasin  drainage  from  the 
Wisconsin  to  the  Fox  through  the  Portage  canal  (a  canal 
15.24  m wide  and  0.61  m deep,  that  connects  the  two  rivers), 
Duck  Creek,  and  low-lying  areas  both  north  and  south  of  Por- 
tage has  occurred  during  high-water  periods.  Levees  along 
the  banks  of  the  Wisconsin  River  near  Portage  protect  that 
city  during  most  flood  flows.  The  levees  also  prevent  any 
widespread  drainage  of  Wisconsin  River  waters  into  the  Fox 
watershed.  The  Columbia  power  plant  could  have  an  effect  on 
the  drainage  situation,  because  of  the  backwater  effect 
created  by  the  power  plant's  encroachment  on  the  Wisconsin 
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River  floodplain. 

The  Columbia  Generating  Station  is  located  on  an  11  km^ 
site  in  the  towns  of  Dekorra  and  Pacific  on  the  east  flood- 
plain  of  the  Wisconsin  River  (See  Figure  1).  The  land  com- 
prising the  site  is  bounded  on  the  south  by  county  trunk 
highway  VJ  and  Rocky  Run  Creek;  on  the  east  by  railroad  tracks 
belonging  to  the  Chicago,  Milwaukee,  St.  Paul,  and  Pacific 
Railroad;  on  the  west  by  the  Wisconsin  River;  and  on  the  north 
by  the  north  section  line  of  Sections  21  and  22,  Town  of  Paci- 
fic, T12N,  R9E  (EIS , 1974). 

The  site  is  in  an  area  of  glacial  drift  overlying  sand- 
stone bedrock.  The  glacial  drift,  an  aquifer,  consists  pri- 
marily of  sand  and  gravel  deposits  with  intermingled  lenses 
of  clay,  silt,  and  coarse  gravel.  Surficial  materials  above 
the  sandstone  also  include  glacial  lake  deposits  averaging 
30.48  m in  thickness.  A series  of  test  wells  in  the  area 
show  depth  to  the  bedrock  of  Upper  Cambrian  sandstone  some- 
where between  21  m and  43  m.  Precambrian  crystalline  rocks 
lie  beneath  the  sandstone.  The  surficial  material  on  the 
floodplain  between  the  Wisconsin  River  and  the  west  dike  of 
the  power  plant's  cooling  lake  is  heterogeneous,  with  sub- 
stantial vertical  and  horizontal  variations.  The  surficial 
glacial  deposts  have  been  modified  by  fluvial  action  and  sur- 
face weathering.  Sediment  deposition  on  the  floodplain  has 
occurred  over  the  western  portion  of  the  site.  Some  windblown 
deposits  exist  in  the  eastern  site  area  (EIS,  1974). 
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An  important  feature,  for  purposes  of  the  present  study, 
is  the  power  plant's  cooling  lake.  The  lake  has  an  approxi- 
mate  surface  area  of  1.94  km'1  and  an  average  water  depth  of 
2.74  m.  The  cooling  lake  normally  contains  about  4.67  bil- 
lion liters  of  water.  Dikes,  composed  of  compacted  fill  of 
local  materials,  surround  the  lake.  Rip-rap  and  bentonite 
cover  the  cooling  lake  sides  of  the  dikes.  A thin  organic 
soil  layer  and  sparse  grasses  top  the  compacted  fill  on  the 
outside  of  the  dikes. 

The  floodplain  surface,  west  of  the  cooling  lake,  and 
the  floodplain  side  of  the  cooling  lake's  west  dike  are  the 
principal  areas  under  investigation  for  the  possibility  of 
scour.  Specifically,  the  area  in  the  vicinity  of  cross  sec- 
tion 8 (See  Figure  5)  is  the  place  of  interest.  Here  a nar- 
row floodplain  channel  exists  between  the  west  dike  and  the 
sandy  knoll  to  the  west  of  the  dike.  Preliminary  Soil  Con- 
servation Service  (SCS)  maps  for  Columbia  County  show  the 
floodplain  surface  in  this  area,  a sedge  meadow,  to  have  a 
Lapeer  fine  sandy  loam  soil.  The  soil  on  the  sandy  knoll, 
labeled  the  North  Knoll  by  the  Columbia  Site  Impact  Study 
Group,  is  a Plainfield  loamy  fine  sand  (EIS,  1974).  The  local 
materials  which  compose  the  compacted  fill  in  the  west  dike 
are  similar  to  the  soil  on  the  North  Knoll,  because  the  fill 
was  extracted  from  this  and  similar  sandy  knolls  nearby.  Sur- 
face soils  currently  found  in  the  sedge  meadow  between  the 
west  dike  and  the  North  Knoll  are  primarily  either  sand  or 
peat,  distributed  in  a non-uniform  fashion  (EIS,  1974). 


Figure  5:  Study  Site  Vegetation  and  Soils  (from  IES,  1976) 
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CHAPTER  III 
GENERAL  PROCEDURE 

Three  hypotheses  are  tested  in  the  present  research. 

The  hypotheses,  expressed  in  the  null  form,  are: 

1)  That  the  increased  flow  velocities,  resulting  from  the 
presence  of  the  Columbia  power  plant,  will  not  cause  scour- 
ing of  the  floodplain  or  west  dike  of  the  cooling  lake  in 
the  vicinity  of  cross-section  8. 

2)  That  the  backwater  effect,  created  by  the  power  plant's 
encroachment  on  the  floodplain,  will  not  produce  an  increase 
in  water  surface  elevation  greater  than  15.24  cm  during  the 
100-year  (regional)  flood. 

3)  That  the  climate  in  the  Upper  Wisconsin  River  Valley  does 
not  have  the  characteristics  of  a predominantly  meridional 
regime. 

If  it  is  determined  that  the  present  climate  resembles 
the  meridional  pre-1895  climatic  regime,  hypotheses  1 and  2 
will  be  re-examined.  Climatic  and  hydrologic  data  which 
should  represent  the  actual  climate  will  be  used  if  possible. 

THE  SCOUR  HYPOTHESIS:  METHOD 
The  effect  of  increased  flow  velocities  may  include 
scour  of  the  floodplain  surface  and  structures.  Scour  can 
occur  if  the  velocity  increase  permits  the  achievement  of  a 
threshold  value  for  some  variable  which  can  initiate  parti- 
cle motion.  The  use  of  the  threshold  variable,  critical  ve- 
locity, is  appropriate  if  one  assumes  that  the  floodplain 
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and  west  dike  surfaces  possess  specific  critical  velocity 
thresholds,  rather  than  intrinsic  geomorphic  thresholds  of 
the  type  suggested  by  Schuiran  (1973). 

The  present  research  employs  a Soil  Conservation  Ser- 
vice method  to  determine  critical  velocity  and  the  possi- 
bility of  scour.  The  method  is  the  "allowable  velocities 
procedure"  for  the  design  of  unlined  earth  channels  (SCS, 
1964,  pp.  6-1  to  6-9).  The  procedure  employs  charts  to  de- 
termine the  non-scouring  velocities  that  can  be  tolerated  by 
the  soil  before  scour  will  occur. 

In  order  to  use  the  method,  certain  soil  and  water  char- 
acteristics must  be  known.  Soil  samples,  taken  from  the 
sedge  meadow  and  west  dike,  are  analyzed  for  soil  classifi- 
cation, plasticity  index,  and  D^,.  (soil  particle  diameter 
greater  than  the  size  of  75%  of  the  sample).  The  soil  clas- 
sification and  are  obtained  by  sieve  and  hydrometer 
analysis.  The  plasticity  index  is  found  by  the  determina- 
tion of  the  water  content,  liquid  limit,  and  plastic  limit 
of  the  samples.  The  plasticity  index  is  simply  the  differ- 
ence between  the  liquid  limit  and  the  plastic  limit  (Lambe, 
1951) . The  character  of  the  water  flowing  across  the  sur- 
faces of  the  floodplain  and  the  west  dike  was  determined  by 
personal  observation  and  grab  sampling  of  suspended  sediment. 
Character  refers  to  what  the  water  carries,  whether  it  be 
suspended  load,  bedload,  or  no  detritus  at  all. 

The  SCS  method  is  designed  for  unvegetated  channels  and 
banks.  The  west  dike,  nearly  covered  with  bluegrass  and 
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sedges,  and  the  floodplain  surface,  a sedge  meadow,  are  vege- 


tated. To  compensate  for  this  fact,  the  assumption  is  made 
that  the  west  dike  and  nearby  floodplain  surface  are  unvege- 
tated. If  the  SCS  method  indicates  the  areas  under  consid- 
eration have  no  potential  for  scour  in  an  unvegetated  condi- 
tion, it  is  apparent  that  they  cannot  scour  in  a vegetated 
state.  Vegetation  strengthens  soil  structure  and  increases 
resistance  to  flow.  Thus,  flow  velocities  and  potential  for 
scour  are  reduced.  If  the  dike  or  floodplain  surface  is  found 
to  be  prone  to  scour,  an  SCS  method  that  considers  the  vege- 
tation can  be  used  to  investigate  the  situation.  However, 
with  this  method,  which  is  based  upon  retardance  factors 
(SCS,  1954),  only  qualitative  estimates  of  the  possibility  of 
scour  can  be  made. 

The  question  of  scour  now  centers  on  the  critical  veloc- 
ity. A somewhat  theoretical  approach  is  used  in  order  to  ar- 
rive at  velocity  values  that  could  be  expected  to  occur  along 
the  dike  and  floodplain  in  the  vicinity  of  cross-section  8. 

The  approach  entails  the  use  of  a Corps  of  Engineers  compu- 
ter program,  the  Hydrologic  Engineering  Center-2  (HEC-2) 
backwater  surface  profile  model  (C  of  E,  1973).  The  program 
can  be  used  to  obtain  velocity  distributions  across  the 
length  of  individual  cross-sections,  to  include  the  veloc- 
ities on  the  floodplain.  The  velocity  values  on  the  flood- 
plain  in  the  reach  between  the  North  Knoll  and  the  west  dike 
near  cross-section  8 can  be  compared  to  the  non-scouring  ve- 
locities on  the  SCS  charts.  If  the  flow  velocities  exceed 
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the  allowable  velocities  determined  by  the  SCS  (1964)  method, 
scour  could  occur. 

Since  the  HEC-2  program  is  so  vital  to  the  present 
study,  a brief  discussion  of  the  program  and  of  its  under- 
lying theories  and  limitations  is  appropriate.  The  program 
computes  and  plots  the  water  surface  profile  for  river  chan- 
nels and  floodplains  for  either  subcritical  or  supercritical 
flow.  The  effects  of  various  hydraulic  structures,  such  as 
bridges,  culverts,  and  dams,  may  be  considered.  Modified 
flow  conditions  resulting  from  channel  improvements,  levees, 
and  floodways  can  also  be  allowed  for. 

Given  the  dimensions  of  each  cross-section,  a designated 
discharge,  and  an  estimated  or  known  value  for  the  starting 
water  surface  elevation  at  the  initial  cross-section,  the 
model  can  calculate  water  surface  elevations  and  flow  veloc- 
ities at  each  cross-section  for  the  channel  and  floodplain. 
The  program  applies  Bernoulli's  Theorem  for  the  total  energy 
and  the  Manning  Equation  for  the  velocity  at  each  cross- 
section  (C  of  E,  1973).  The  model  applies  these  relation- 
ships to  successive  cross-sections  until  the  entire  reach 
under  investigation  is  accounted  for.  The  essential  equa- 
tions are  here  presented  in  order  to  clarify  the  discussion. 
The  Manning  Equation: 

V = 1 R.2/ 3 s1/2, 
n 

where  V,  the  velocity  (unknown),  equals  one  divided  by  n, 
the  Manning  roughness  coefficient  (estimated),  times  r2/3> 
the  hydraulic  radius  (known)  to  the  2/3  power,  times  sV2f 
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the  riverbed  slope  (known)  to  the  1/2  power. 

Bernoulli's  Equation: 

Vl2  + Dl  + = V22  + D2  + Z2  + hL 

2g  2g 

expresses  the  energy  relationships  between  two  cross-sections 

(1  and  2),  where  V^2,  the  velocity  head  (V^2  is  velocity  from 

2g 

the  Manning  Equation  squared;  it  is  divided  by  twice  the  gravi- 
tational constant)  at  cross-section  1,  plus  Dp  the  hydraulic 
gradient  or  depth  at  cross-section  1 (known),  plus  Zp  the 
riverbed  elevation  at  cross-section  1 (known),  equals  the  sum 

of  the  values  for  the  same  variables  for  cross-section  2 
2 

(\?2  /2g  + D2  + Z2)  plus  hp  the  head  loss.  Generally,  the 

head  loss  is  due  to  friction,  so  h^  = hp  (friction  head  loss). 
The  computer  model  computes  the  friction  head  loss  with  the 
relationship : 
hp  = S2  + A L, 

where  S2  and  are  the  slopes  of  the  energy  gradient  (ob- 
tainable from  Bernoulli's  Equation,  since  S = V2/2g  + D + Z) 
at  the  two  cross-sections  and  AL  is  the  length  of  the  reach 
between  the  two  cross-sections  (known). 

From  Bernoulli's  energy  relationship,  it  is  apparent  that 
all  the  terms  on  either  side  of  the  equation,  except  for  D2 
(D^  is  known  since  it  is  the  estimated  or  otherwise  deter- 
mined starting  elevation  for  the  initial  cross-section) , are 
known.  The  equation  can  be  solved  for  the  water  depth  at 
cross-section  2 (D2) , because  D2  is  the  only  unknown.  Thus, 
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the  Manning  Equation  and  Bernoulli's  Theorem  enable  the 
model  to  calculate  velocity  and  water  surface  elevation, 
which  are  the  variables  needed  in  a study  of  scour  and  back- 
water problems.  The  following  diagram  depicts  the  energy 
relationship  expressed  by  Bernoulli's  Equation. 


Figure  6:  Bernoulli  Energy  Relationships 


Horizontal  line 


(from  Simons,  1969,  p.  134) 


The  HEC-2  model  is  widely  employed  by  various  agencies 
throughout  the  country,  including  the  Corps  of  Engineers, 
Wisconsin  DNR,  and  private  engineering  firms.  The  model  is 
recognized  as  an  outstanding  tool  when  used  with  accurate 
data  and  interpreted  with  objectivity  (Hampton,  1976). 

Other  agencies,  such  as  the  United  States  Geological  Sur- 
vey and  Bureau  of  Public  Roads,  have  their  own  programs. 

' 
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However,  all  programs  are  basically  the  same.  The  only  dif- 
ference between  the  USGS  model  and  the  HEC-2  program  in- 
volves the  computations  of  flow  through  bridges  (Krug,  1976; 
Hampton,  1976).  Since  bridges  do  not  affect  the  flow  in 
the  study  area,  it  is  reasonable  to  assume  that  any  back- 
water profile  model  should  produce  similar  results  if  used 
in  the  present  research. 

Limitations  necessarily  accompany  the  HEC-2  model.  A 
major  limitation  is  cost.  Surveyed  cross  sections  are  vital 
for  accurate  computations.  The  present  reserach  is  fortu- 
nate in  that  sufficient  cross  sections  were  surveyed  for  the 
engineering  and  environmental  studies  prior  to  the  construc- 
tion of  the  power  plant.  Twenty-one  cross  sections,  sur- 
veyed in  the  early  1970s  are  contained  within  the  study 
reach  (See  Figure  2) . The  sections  are  accurately  des- 
cribed with  20  to  80  horizontal  ground  points  and  eleva- 
tions per  section.  Cross-section  data  were  reviewed  for 
accuracy  and  simplification.  Redundant  points  were  elima- 
ted. 

Computer  time  and  operator  salaries  add  to  the  cost 
limitation;  however,  the  wide  applicability  of  the  computer 
output  should  provide  adequate  compensation.  A bigger  limi- 
tation may  be  the  insensitivity  of  the  Manning  "n"  variable, 
the  roughness  coefficient.  The  program  can  account  for 
changes  in  flow  resistance  with  increasing  water  discharge 


and  depth  in  the  river  channel,  but  it  cannot  perform  simi- 
larly for  flow  across  the  floodplain.  On  the  other  hand, 
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the  model  allows  for  multiple  Manning  "n"  values  along  the 
length  of  the  cross-section:  for  example,  cross-section  8 
has  five  different  roughness  coefficients.  Manning  "n" 
values  from  Lee's  original  DNR  work  on  the  backwater  prob- 
lem at  Portage  were  used  in  the  present  analysis.  The 
roughness  values  were  checked  against  air  photographs  and 
on-site  observations  for  accuracy,  but  for  the  most  part, 
reliance  on  the  judgement  of  Mr.  Lee  and  his  co-workers  was 
necessary. 

Another  limitation  of  the  model  is  the  need  to  estab- 
lish a starting  elevation  at  the  initial  cross-section. 

For  the  100-year  flood,  the  starting  elevation  should  be 
quite  accurate,  since  it  was  obtained  from  the  Corps  of  En- 
gineers' floodplain  information  study  for  Portage,  Wiscon- 
sin (C  of  E,  1972).  However,  the  starting  elevations  for 
floods  of  other  recurrence  intervals  were  dependent  upon 
the  engineering  judgement  of  Mr.  Terry  Hampton  of  the  DNR 
Floodplain-Shoreland  Management  Section.  Mr.  Hampton  ran 
the  HEC-2  computer  program  for  the  present  study. 

The  only  other  input  needed  for  the  HEC-2  model  is  the 
water  discharge.  Discharge  values  from  Conger  (1971),  as 
modified  by  the  USGS  for  the  2,  10,  and  100-year  floods  at 
Wisconsin  Dells,  were  used.  There  is  no  gaging  station  at 
Portage,  only  a crest  gage,  necessitating  the  use  of  Wiscon- 
sin Dells  data.  The  program  was  run  for  all  three  floods 
with  cross-sections  that  existed  before  the  Columbia  power 
plant  was  constructed.  The  outputs  from  the  runs  are  the 
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"natural  state"  data  which  serve  as  the  bases  for  compari- 
son with  output  from  runs  which  consider  the  power  plant's 
encroachment.  Program  runs  were  made  for  the  present 
situation,  that  is,  with  cross-sections  affected  by  the 
power  plant's  dikes,  and  for  the  situation  which  considers 
double  encroachment.  In  the  double  encroachment  runs,  the 
amount  of  conveyance  lost  because  of  the  presence  of  the 
plant  is  subtracted  from  the  conveyance  on  the  opposite 
floodplain  if  the  physical  situation  allows  it.  The  same 
sequence  of  runs  was  followed  for  all  three  recurrence  in- 
tervals. The  velocity  output  for  cross-section  8 provides 
the  mean  velocity  values  which  should  occur  along  the  flood- 
plain  surface  and  west  dike  for  each  of  the  three  floods 
under  consideration.  A comparison  of  these  values  with  the 
allowable  velocities  for  the  floodplain  and  west  dike  soils, 
as  determined  by  the  SCS  (1964)  method,  allows  one  to  test 
the  first  hypothesis-- that  increased  flow  velocities  re- 
sulting from  the  power  plant's  encroachment  will  not  cause 
scour  of  the  floodplain  or  west  dike  in  the  vicinity  of 
cross-section  8. 

THE  BACKWATER  HYPOTHESIS:  METHOD 
Analysis  of  output  from  the  HEC-2  computer  model  de- 
termines the  increases  in  water  surface  elevation  because  of 
the  backwater  effect  caused  by  the  power  plant's  encroach- 
ment. The  program  computes  the  water  surface  elevation  for 
each  cross  section.  Elevation  differences  between  the 
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natural  state  runs  and  the  runs  for  double  encroachment  are 
results  of  the  proposed  backwater  effect.  The  elevation 


differences  at  each  cross  section  are  the  increases  in 
water  surface  elevation  resulting  from  encroachment. 

The  backwater  effect  is  of  significance,  not  only  be- 
cause the  increase  in  water  surface  elevation  may  be  large 
enough  to  be  in  violation  of  the  state  floodplain  zoning 
ordinance,  but  also  because  of  its  possible  impact  on  the 
city  of  Portage.  Some  flood  protection  is  afforded  the  city 
by  the  levees  along  the  banks  of  the  Wisconsin  River;  how- 
ever, the  levees  are  only  expected  to  withstand  floods  with 
recurrence  intervals  of  less  than  10  years  (C  of  E,  1971). 
Larger  floods  with  stages  greater  than  242.32  m at  Portage 
will  overtop  the  levees.  The  100-year  flood,  with  a dis- 
charge of  2690.4  m^/s,  is  expected  to  have  a river  stage  of 
242.93  m.  The  largest  flood  on  record  at  Portage,  a Sep- 
tember 1938  flood  with  a discharge  of  approximately  2044.7 
m^/s  and  a recurrence  interval  of  about  10  years,  reached 
a stage  of  242.50  m.  Fortunately,  the  levees  were  not  over- 
topped because  an  early  warning  allowed  the  citizens  to 
heighten  them  with  sandbags  (C  of  E,  1971).  None  of  the 
aforementioned  figures  consider  the  backwater  effect  created 
by  the  power  plant. 

A further  complication  for  inhabitants  of  Portage  is 
the  fact  that  the  existing  levees  do  not  meet  current  Corps 
of  Engineers  and  DNR  standards  (C  of  E,  1971).  The  Corps 
of  Engineers  estimates  that  annual  damages  from  hypothetical 
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floods  exceeding  levee  design  height  would  average  $225,000 
for  Portage  (C  of  E,  1971).  Preliminary  improvement  plans 
have  been  formulated,  and  some  money  is  presently  available 
for  initial  work  on  flood  protection.  The  Corps  of  Engi- 
neers, DNR,  and  the  Department  of  Housing  and  Urban  Develop- 
ment (HUD)  are  involved,  but  the  end-product  of  their  ef- 
forts appears  to  be  far  off  in  the  future. 

It  is  quite  possible  that  the  encroachment  of  the  power 
plant  increases  the  water  surface  height  to  such  an  extent 
that  floods  more  frequent  than  the  10-year  event  might  over- 
top the  levees.  An  analysis  of  the  HEC-2  output  will  not 
only  test  the  present  study's  second  hypothesis,  that  the 
backwater  effect  will  not  create  a greater  than  15.24  cm 
water  surface  elevation  increase  at  Portage  during  the  100- 
year  flood,  but  it  will  also  evaluate  the  significant  prob- 
lem posed  above. 

THE  CLIMATIC  CHANGE  HYPOTHESIS:  METHOD 

A global  cooling  with  accompanying  rainfall  anomalies 
since  about  1960  is  widely  accepted  as  fact.  Lamb's  hypothe- 
sis, that  the  cooler  and  wetter  regime  of  today  is  a result 
of  a shift  in  the  tropospheric  long  wave  patterns  (Lamb, 

1966) , is  a reasonable  statement  for  the  recent  climatic 
change.  Knox  and  other  (1975)  investigated  the  situation 
in  the  Upper  Mississippi  River  Valley  and  found  the  pro- 
posed change  to  be  quite  evident. 


The  present  study  analyzes  precipitation  and  stream- 
flow  data  for  the  Upper  Wisconsin  River  Valley  from  the  Lac 
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Vieux  Desert  to  the  Columbia  power  plant  in  an  attempt  to 
test  for  such  a change  there.  Rainfall  and  low  flow  stages 
for  the  period  1873-1910  and  precipitation  and  mean  dis- 
charge per  square  kilometer  for  the  interval  1925-1975  are 
the  data  for  the  investigation.  The  data  are  evaluated  in 
order  to  determine  the  presence  or  absence  of  a trend  to- 
ward a greater  degree  of  wetness  and  of  a persistence  of 
variability  in  rainfall  and  hydrologic  response  for  the  pre- 
1895  and  post-1950  periods.  Persistent  variability  results 
from  the  alternating  flow  of  air  from  the  north  and  the 
south,  as  the  meridional  circulation  pattern  blocks  the 
westerly  airflow  of  the  zonal  regime.  A trend  toward  in- 
creased precipitation  would  result  primarily  from  the  south- 
erly component  of  the  meridional  regime.  The  incursion  of 
warm,  moist  air  from  the  Gulf  of  Mexico  enhances  precipita- 
tion in  regions  like  the  Upper  Wisconsin  River  Valley. 

The  study  of  climate  change  is  divided  into  two  parts. 
The  first  pertains  to  the  1873-1910  data,  while  the  second 
part  is  concerned  with  the  more  recent  data.  Precipitation 
records  for  the  initial  period  were  obtained  from  an  early 
work  by  Mead  (1911),  who  examined  rainfall-runoff  relation- 
ships in  Wisconsin  at  the  turn  of  the  present  century.  He 
used  five  weather  stations  located  upstream  of  Portage  to 
obtain  his  monthly  rainfall  totals  for  the  Upper  Wisconsin 
River  Valley.  Monthly  low  flows  at  Portage  are  the  stream- 
flow  data  for  the  first  part  of  the  climate  study.  Low 
flows  are  appropriate,  because  the  Wisconsin  River  was  not 


regulated  to  any  great  extent  during  the  1873-1910  period. 
The  low  flows  should  be  good  indicators  of  hydrologic  re- 
sponse to  precipitation  and,  thus,  to  the  climatic  regime. 


The  rainfall  and  streamflow  data  sets  are  grouped  into 
three  hydrologic  seasons.  Monthly  groupings  of  November  to 
February,  March  to  June,  and  July  to  October  represent  the 
Winter,  Spring,  and  Fall  seasons  respectively  (Knox  et  al . , 
1975).  The  rainfall  and  streamflow  data  are  also  expressed 
as  running  means  and  cumulative  deviations  from  the  mean  in 
order  to  detect  any  trends  or  persistent  variability  in  the 
data. 

The  second  part  of  the  investigation  of  climate  change 
involves  the  analysis  of  more  recent  data.  Precipitation 
totals  from  13  weather  stations  (Table  1)  are  examined  for 
the  1925-1975  period.  As  in  the  first  part  of  the  study, 
tve  data  are  divided  into  the  three  hydrologic  seasons. 

The  monthly  precipitation  totals  for  the  13  stations  are 
averaged  to  provide  one  rainfall  value  per  month  for  the  en- 
tire area  above  the  Columbia  power  plant.  The  average  rain- 
fall values  for  each  hydrologic  season  are  presented  in  a 
time  series.  Running  means  and  cumulative  deviations  from 
the  mean  are  employed  for  the  same  purposes  as  in  the  first 
part  of  the  analysis  of  climate  change. 

Streamflow  data  differ  from  those  used  for  the  1873- 
1910  period.  The  Wisconsin  Valley  Improvement  Authority 
received  a charter  in  1907  to  regulate  the  flow  of  the  Wis- 
consin River  with  a series  of  reservoirs.  Power  dams  and 
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Table  1:  Weather  Stations  and  Periods  of  Record 

1.  Phelps  Deerskin  Dam  (1945-  ) 

2.  Rest  Lake  (1913-  ) 

3.  Long  Lake  Dam  (1908-  ) 

4.  Minoqua  Dam  (1903-  ) 

5.  Rhinelander  (1891-  ) 

6.  Tomahawk  Spirit  Reservoir  (1902-  ) 

7.  Merrill  (1905-  ) 

8.  Marshfield  Experimental  Farm  (1912-  ) 

9.  Stevens  Point  (1893-  ) 

10.  Wisconsin  Rapids  (1921-  ) 

11.  Coddington  (1921-  ) 

12.  Hancock  Experimental  Farm  (1902-  ) 

13.  Mauston  (1891-  ) 


Table  2:  Gaging  Stations,  Areas,  and  Periods  of  Record 

1.  Wisconsin  R.iver  at  Rainbow  Lake  1942.5  km^  (1936-  ) 

2.  Tomahawk  River  at  Bradley  1411.6  km2  (1930-September  1973) 

3.  Spirit  River  at  Spirit  Falls  212.4  km2  (1942-  ) 

4.  Prairie  River  near  Merrill  468. km^  (1939-  ) 

5.  Wisconsin  River  at  Merrill  7200.2  km^  (1902-  ) 

6.  Eau  Claire  River  at  Kelly  844.3  km2  (1914-26,  1939-  ) 

7.  Wisconsin  River  at  Rothschild  10360  km2  (1944-  ) 

8.  Big  Eau  Pleine  River  near  Stratford  580.2  km2  (1914-25, 

37-  > 

9.  Yellow  River  at  Babcock  577.6  km2  (1944-  ) 

10.  Lemonweir  River  at  New  Lisbon  1295  km2  (1944-  ) 

11.  Wisconsin  River  at  Wisconsin  Dells  20279.7  km2  (1934-  ) 

12.  Baraboo  River  near  Baraboo  1554  km2  (1913-22,  1942-  ) 

Total  Area:  46,725.5  km3 
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additional  reservoirs  added  to  the  degree  of  regulation  in 
later  years.  Lew  Sherrar  (1976b),  a manager  of  the  regula- 
tory agency,  recommended  the  use  of  monthly  mean  flow  per 
unit  area  for  the  1925-1975  period.  Mean  discharges  are 
used  because  river  regulation  provides  unnatural  control  on 
the  low  flows.  Mean  monthly  discharges  for  12  gaging  sta- 
tions (Table  2)  in  the  Upper  Wisconsin  River  Valley  above 
Portage  are  examined  for  the  1941-1973  interval.  The 
stream  record  is  not  long  enough  to  employ  the  same  1925- 
1975  interval  used  for  the  precipitation  analysis.  The 
streamflow  data  are  presented  and  analyzed  in  the  same 
manner  as  are  the  precipitation  data. 

A final  measure  used  to  examine  the  possibility  of  a 
change  in  climate  is  an  analysis  of  the  temporal  and  abso- 
lute differences  in  the  averages  of  monthly  rainfall  for 
each  of  13  weather  stations.  The  monthly  rainfall  data  for 
each  station  are  divided  into  a 1925-1949  zonal  interval 
and  a 1950-1975  meridional  interval  for  the  purpose  of  com- 
parison. 
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CHAPTER  IV 

THE  SCOUR  HYPOTHESIS:  DATA  AND  ANALYSIS 

Basically,  the  scour  investigation  considers  two  vari- 
ables: stream  velocity  and  soil  character.  Stream  veloci- 
ties are  determined  with  the  HEC-2  backwater  surface  pro- 
file computer  model.  Various  soil  analysis  techniques  pro- 
vide data  on  the  soil  characteristics. 

Velocities 

The  levees  at  Portage  present  a complication  in  the 
HEC-2  analysis.  The  Corps  of  Engineers  expects  the  levees 
to  fail  during  high-magnitude  flood  events,  whereas  Mr. 
Hampton  of  the  DNR  suspects  that  they  may  not  (C  of  E,  1971; 
Hampton,  1976) . In  the  upper  reach  of  the  study  area  near 
Portage,  the  latter  situation  would  confine  most  of  the 
floodwater  inside  the  levees.  The  present  research  con- 
siders both  hydraulic  settings  in  the  scour  analysis  (Fig- 
ures 7A  and  B) . 

The  water  flow  concept  of  the  Corps  of  Engineers  for 
extreme  floods  is  depicted  in  Figure  7A.  The  levees  will 
fail,  allowing  water  to  flow  across  the  length  of  each 
cross-section  if  the  volume  of  water  is  sufficient  to  en- 
compass the  area  at  each  section.  The  DNR's  concept  (Fig- 
ure 7B)  confines  the  floodwater  to  the  levees  at  cross- 
sections  16,  14,  and  13.  Water  at  stages  greater  than  levee 
height  overflows  and  possibly  moves  into  storage  or  drains 
east  to  the  Fox  River  watershed.  Because  of  the  confine- 


Figure  7A:  Corps  of  Engineers  Flow  Concept 
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ment  of  flow  near  Portage,  the  DNR's  hydraulic  setting  is 
the  more  streamlined  in  the  middle  reaches.  Flow  will  not 
occur  across  the  total  lengths  of  each  cross-section,  be- 
cause of  the  cut-offs  at  the  levees  and  the  streamlined 
flow  situation  below  the  levees. 

Three  flood  discharges  are  examined  under  three  dif- 
ferent conditions:  the  "natural  state"  that  existed  before 
the  construction  of  the  power  plant,  the  present  situation 
with  the  Columbia  plant  encroachment,  and  the  legal  condi- 
tion which  considers  double  encroachment  (the  requirement  to 
consider  the  effects  of  an  identical  project  on  the  other 
side  of  the  river) . A comparison  of  computer  output  for 
the  first  two  conditions  enables  one  to  determine  the  ve- 
locity increases  and  the  decreases  in  the  amount  of  dis- 
charge conveyed  on  the  east  floodplain  as  a result  of  the 
power  plant's  encroachment  (Table  3). 

The  velocity  values  for  the  present  condition  at  cross- 
section  8 are  the  data  needed  to  investigate  the  scour  po- 
tential. The  HEC-2  output  contains  the  mean  velocities  ex- 
pected to  exist  in  the  sedge  meadow  and  along  the  face  of 
the  west  dike  (Table  4).  Cross-section  8 contains  41  hori- 
zontal ground  points  and  elevations.  The  first  ground 
point,  at  horizontal  station  zero,  is  the  left  end  point  of 
the  cross-section  as  one  looks  downstream.  The  top  of  the 
west  dike  is  located  at  the  487.68  m horizontal  ground 
point,  and  the  east  side  of  the  North  Knoll  is  located  near 
the  690.60  m ground  point.  Velocity  values  that  exist  be- 


Table  3:  Velocity  Increases  and  Conveyed  Discharge  Decreases  Near  West  Dike  Attributable 
to  Encroachment  of  Columbia  Power  Plant. 
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Table  4:  Velocities  and  Percent  Total  Discharge  Conveyed 
Within  Specific  Ground  Points  at  Cross-sections 
Near  Columbia  Power  Plant. 

DNR  SETTING 

Q = 979.9  m^/s  (2-year  flood) 

Cross-  Ground-  Natural  Present 


sections 

points  (in) 

v(m/s) 

7oQ 

v(m/s) 

7oQ 

8 

142.0-  487.7 

.15 

4.9 

487.7-  666.9 

.18 

3.9 

.18 

4.1 

666.9-  777.5 

.18 

3.1 

.21 

3.3 

7 

304.8-  670.6 

.12 

3.5 

__ 

670.6-  914.4 

.15 

4.5 

_ — 

_ _ 

829.4-  975.4 

-- 

-- 

.15 

3.2 

6.5 

2194.6-2575.6 

.09 

5.7 

.09 

5.7 

Q = 1761.5  m^/s  (10-year  flood) 


138.1- 

304.8 

.21 

3.8 

— Mi 

_ _ 

304.8- 

487.7 

.24 

5.4 

— 

— — 

487.7- 

609.6 

.27 

3.9 

.27 

4.3 

609.6- 

731.5 

.27 

4.5 

-.30 

4.9 

304.8- 

487.7 

.15 

3.0 

__  — 

487.7- 

670.6 

.18 

3.7 

— CM 

— — 

670.6- 

914.4 

.18 

6.2 

- — 

_ Ml 

829.4- 

975.4 

.21 

4.4 

2194.6- 

2575.6 

. 12 

7.5 

.12 

7.5 

Q = 2690.4 

3. 

m /s 

(100-year 

flood) 

8 

134.4-  243.8 

.30 

3.1 

__ 

243.8-  487.7 

.30 

8.7 

- — 

— — 

487.7-  609.6 

.34 

4.7 

.37 

5.3 

609.6-  731.5 

.34 

5.1 

.40 

5.8 

7 

304.8-  452.3 

.21 

3.1 

__ 

452.3-  609.6 

.21 

3.7 

— — 

— _ 

609.6-  914.4 

.24 

8.5 

- — 

_ _ 

829.4-  975.4 

-- 

.27 

5.0 

6.5 

2194.6-2575.  6 

.15 

8.5 

.15 

8.5 

Table  4:  continued 

CORPS  OF  ENGINEERS  SETTING 

Q = 979.9  m^/s  (2-year  flood) 

Cross-  Ground-  Natural  Present 

sections points  (m) V(m/s)  7»Q V(m/s)  7<,Q 
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8 

142.0-  487.7 

.15 

5.0 

- — 

487.7-  666.9 

.18 

3.9 

.18 

4.1 

666.9-  777.5 

.18 

3.1 

.21 

3.3 

7 

4.3-  548.6 

.09 

3.3 

_ _ 

548.6-  914.4 

.15 

6.0 

— 

— 

829.4-  975.4 

-- 

-- 

.15 

3.2 

6.5 

555.3-1447.8 

.06 

5.3 

_ — 

' 

1045.5-1447.8 

- - 

- “ 

.06 

3.8 

Q = 1761.5 

m^/s 

(10  -year 

flood) 

8 

138.1-  304.8 

.21 

3.8 

_ — 

-- 

304.8-  487.7 

.24 

5.4 

- - 

- 

487.7-  609.6 

.27 

3.9 

.27 

4.3 

609.6-  731.5 

.27 

4.5- 

.30 

4.9 

7 

(-) 30. 8-  304.8 

.15 

3.8 

_ _ 

304.8-  548.6 

.15 

4.0 

- - 

— 

548.6-  914.4 

.18 

8.4 

- - 

— 

829.4-  975.4 

-- 

-- 

.21 

4.4 

6.5 

320. 6-1045.5 

.06 

4.1 

_ — 

_ — 

1045.5-1447.8 

.09 

5.5 

.09 

5.7 

Q = 2690.4 

m^/s 

(100-year 

flood) 

8 

134.4-  243.8 

.30 

3.0 

_ _ 

_ 

234.8-  487.7 

.30 

8.6 

- - 

- - 

487.7-  609.6 

.34 

4.7 

.37 

5.3 

609.6-  731.5 

.34 

5.1 

.40 

5.8 

7 

(-) 86 . 0-  182.9 

.15 

3.1 

_ _ 

_ _ 

182.9-  365.8 

.21 

3.6 

- - 

— 

365.8-  548.6 

.21 

3.8 

-- 

548.6-  914.4 

.21 

9.4 

— 

— 

829.4-  975.4 

-- 

-- 

.27 

5.0 

6.5 

239.0-1045.5 

.09 

6.4 

_ _ 

_ _ 

L_ 

1045.5-1447.8 

.12 

6.3 

.12 

6.7 
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tween  these  last  two  points  are  the  mean  velocities  to  be 
compared  with  the  allowable  velocities  obtained  with  the 
SCS  critical  velocity  method. 

Soils 

Soil  samples  from  two  distinct  areas  are  analyzed  in 
the  present  study.  It  must  be  emphasized  that  the  soils  are 
considered  in  the  engineering  sense,  that  is,  any  material 
from  the  surface  down  to  bedrock  is  soil.  The  floodplain 
surface  (a  sedge  meadow)  near  cross-section  8 contains  peat, 
marl,  and  sand.  The  dominant  surface  soil  in  the  sedge 
meadow  is  peat,  averaging  nearly  60  cm  in  depth  and  over- 
lying  marl  and  fine  sand.  Peat  depths  increase  up  to  3.05  m 
near  cross-sections  7 and  6.5  and  up  to  1.52  m between 
cross-sections  8 and  8.1  (IES,  1976). 

Medium  and  fine  sands  are  exposed  on  the  floodplain 
surface  in  two  specific  locations:  the  bed  of  a drainage 
channel  which  flows  north  through  the  sedge  meadow  to  Duck 
Creek,  and  the  center  of  the  study  area  where  an  access  road 
once  existed  and  the  sedge  vegetation  is  sparse  (Figure  5) . 
The  latter  area  contains  some  coarse  sand  and  gravel  (pos- 
sibly remnants  of  the  access  road  and  associated  fill  ma- 
terials) . The  sedge  meadow  surface  soils  are  obviously  not 
members  of  the  same  population;  therefore,  only  five  sam- 
ples, thought  to  be  representative  of  the  peat  and  sand  de- 
posits, are  analyzed  (Table  5). 


A 


Table  5:  Soil  Sample  Characteristics 
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# dike  subgroup  near  cross-section 


The  west  dike  is  composed  of  compacted  fill  throughout 
its  structure  (EIR,  1973).  Thirty  soil  samples,  obtained 
in  a systematic  fashion  from  the  dike  (Figure  8) , are  analy- 
zed for  their  classification,  D75,  and  plasticity.  Since 
the  dike  is  composed  of  a homogeneous  distribution  of  soil, 
representative  data  for  the  description  of  the  west  dike  are 
obtainable . 

The  30  samples  from  the  dike  and  the  5 samples  from  the 
floodplain  are  surface  samples,  generally  440  cm3  in  s±Ze. 
Entire  samples,  air-dried,  were  used  in  the  sieve  analysis 
(Appendix  1)  in  order  to  eliminate  bias.  Only  50  grams 
(smaller  than  0.5  mm)  of  each  sample  were  used  in  the  hydro- 
meter analysis  (Appendix  2) , and  approximately  30  grams 
(smaller  than  2.0  mm)  were  used  in  the  plasticity  determina- 
tions (Appendix  3).  Cumulative  frequency  curves  for  each 
sample  provide  a graphical  presentation  of  the  values 
(at  25%  on  the  curves)  and  the  percentage  of  sand  and 
coarser  materials  (at  4 phi  on  the  curves)  for  each  sample 
(Figures  9A  and  9B) . Plasticity  index  values  are  obtained 
by  the  determination  of  the  liquid  and  plastic  limits  with 
the  liquid  limit  device  or  the  expedient  method. 

The  floodplain  soils  are  difficult  to  define  because  of 
the  peat  and  marl  deposits.  The  peat  appears  cohesive,  as 
the  deposits  are  either  below,  at,  or  slightly  above  the 
water  taole.  Yet  the  peat  has  no  measurable  plasticity. 

The  liquid  limit  of  the  peat  is  greater  than  1007o,  that  is, 
the  weight  of  water  held  by  a sample  is  greater  than  the 
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weight  of  the  peat  sample  itself.  It  appears  that  the  peat 
deposits  present  a dilemma.  According  to  Brady  (1974), 
peat  has  low  plasticity  and  cohesion,  but  does  not  behave 
as  a mineral  soil  with  the  same  characteristics  because  of 
the  peat's  fibrous  structure.  Peat  normally  has  water  con- 
tents much  greater  than  100%,  in  fact,  peat  can  hold  water 
to  the  extent  of  one  to  20  times  its  weight.  If  the  peat 
dries  out,  it  can  harden  to  a pan-like  structure  or  break 
apart  into  small  particles  subject  to  wind  erosion  (Brady, 
1974).  The  study  area  should  not  reach  such  a state  because 
the  ground  water  level  is  at  the  surface  throughout  most  of 
the  sedge  meadow.  This  condition  has  been,  and  should  con- 
tinue to  be,  aided  by  the  inflow  of  water  from  the  cooling 
lake  either  through  or  under  the  base  of  the  dike  into  the 
sedge  meadow  (IES,  1976).  Thus,  the  peat  deposits  in  the 
sedge  meadow  should  always  be  wet.  In  such  a state  the  peat 
seems  to  offer  a high  resistance  to  scour,  even  though  the 
material  has  low  plasticity  and  cohesion.  Since  the  criti- 
cal velocity  method  depends  upon  particle  size  (among  other 
things) , the  scour  potential  of  the  peat  cannot  be  evaluated 
in  this  research  because  of  the  peat's  fibrous  nature. 
Therefore,  only  the  sand  deposits  in  the  drainage  channel 
and  in  the  center  strip  near  cross-section  8 are  considered. 

The  sand  deposits  on  the  floodplain  can  be  evaluated 
for  their  potential  to  scour.  The  four  sand  samples  pre- 
sent reasonably  uniform  characteristics.  They  have  average 
sand-silt-clay  percentages  of  88-2-10,  which  classifies 
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these  surface  soils  as  Sand  in  the  United  States  Department 
of  Agriculture  (USDA)  scheme.  Dy^  values  range  from  1 phi 
to  2 phi,  with  1 phi  (0.5  mm)  selected  as  the  most  repre- 
sentative size.  These  characteristics  only  describe  the 
sand  deposits  in  the  drainage  channel  and  the  center  strip. 
The  greater  area  of  the  sedge  meadow  contains  the  rather 
elusive--from  the  viewpoint  of  critical  velocity--peat  de- 
posits at  the  surface. 

The  dike  presents  a more  definable  situation.  The  20 
samples  taken  near  cross-section  8 have  average  sand-silt- 
clay  percentages  of  89-6-5  and  a USDA  classification  of 
Sand.  Two  of  the  samples  have  Dy^  values  greater  than  2 
phi,  but  the  remaining  18  samples  have  Dy^  sizes  between  0 
phi  and  2 phi.  Again,  1 phi  (0.5  mm)  is  the  most  represen- 
tative Dy^  value.  The  distribution  of  soil  samples  for  the 
entire  dike,  including  samples  21-25  along  cross-section  7 
and  samples  26-30  from  cross-section  6.5,  has  average  sand- 
silt-clay  percentages  of  84.5-8.5-7.0  and  a USDA  classifi- 
cation of  Sand.  Dy^  values  include  two  3 phi  and  eight  1 
phi  sizes  in  the  additional  10  samples.  A representative 
size  of  1 phi  is  selected  as  the  Dy^  value  for  the  entire 
west  dike. 

The  plasticity  of  the  dike  soils  is  very  low.  Only 
samples  near  the  base  of  the  dike  and  samples  from  cross- 
sections  7 and  6.5  exhibit  any  measurable  plastic  limit 
(Figure  10).  Liquid  limits  were  determined  for  all  dike  and 
floodplain  samples  with  the  liquid  limit  device.  The  field 
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10:  Distribution  of  Plasticity  Indices 
s on  south  end  of  dike  in  fine  print) 
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expedient  method  was  also  used  to  determine  liquid  limits 
for  the  20  dike  samples  from  the  vicinity  of  cross-section 
8. 

Every  soil  has  a definite  value  of  the  liquid  limit. 

This  characteristic  adds  to  the  strength  of  any  soil  clas- 
sification scheme.  Sands  normally  have  liquid  limits  that 
range  between  15 % and  20%.  If  the  liquid  limit  is  less  than 
35%,  the  soil  has  low  plasticity.  Soils  with  low  liquid 
limits  are  softened  when  wetted  and  their  strength  is  often 
reduced  by  only  minor  increases  in  water  content.  However, 
soil  structure  has  to  be  destroyed  before  a soil  is  weak- 
ened enough  to  flow,  even  if  the  water  content  is  increased 
up  to  the  liquid  limit  (Kezdi,  1974). 

The  20  dike  samples  taken  near  cross-section  8 exhibit 
liquid  limits  ranging  from  14.7%  to  23.6%,  with  the  liquid 
limit  device.  The  range  is  from  16.4%  to  30.1%  with  the 
expedient  method.  All  20  samples  have  liquid  limit  values 
less  than  35%,  for  either  method.  The  surface  soils  on  the 
dike  near  cross-section  8 exhibit  low  plasticity,  if  any. 

Plastic  limits  could  not  be  measured  for  any  of  the  20  sam- 
ples; thus,  they  all  have  plasticity  indices  of  zero.  The 
plasticity  index  is  the  difference  between  the  liquid  limit 
and  the  plastic  limit,  but  soils  with  no  plastic  limit  have 
no  plasticity  index  either  (Kezdi,  1974). 

The  samples  taken  from  the  dike  near  cross-sections  7 
and  6.5  have  some  plasticity.  In  fact,  6 of  10  samples  had 
measurable  plastic  limits.  The  plastic  limits  ranged  from 
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15.2%  to  20.8%,  which  is  reasonable  since  very  fine  sands 
exhibit  some  plasticity--usually  values  of  17%  to  207„ 

(Kezdi,  1974).  The  plasticity  indices  for  5 of  the  6 sam- 
ples with  plastic  limits  are  greater  than  7.  Bureau  of  Re- 
clamation results  show  a value  of  7 to  be  a tentative  criti- 
cal value  for  canal  design.  Scour  can  result  from  moderate 
threshold  forces  in  materials  with  a plasticity  index  below 
7 (Chow,  1959). 

The  floodplain  samples,  excluding  the  peat  sample,  are 
rather  homogeneous  with  average  sand-silt-clay  percentages 
of  88-2-10.  They  can  be  classified  as  Sand  according  to  the 
USDA  method.  The  liquid  limits  are  from  14.1%  to  21.2%, 
within  the  range  expected  for  fine  sands. 

Analysis 

The  SCS  critical  velocity  charts  can  be  examined  now 
that  the  expected  mean  velocities  and  the  soil  character- 
istics are  known.  All  samples  from  the  floodplain  (peat  ex- 
cluded) and  the  dike  are  categorized  as  either  silty  or 
clayey  sands  with  more  than  5%  fines  (Figure  11A) . Values 
for  plasticity  index  (PI)  exist  only  for  6 samples  from  the 
southern  portion  of  the  dike  near  cross-sections  7 and  6.5. 
All  other  samples  have  plasticity  indices  of  zero.  The  rep- 
resentative value  of  the  is  0.5  mm  (1  phi)  for  all  sam- 
ples. Grab  sampling  of  suspended  sediment  during  an  ap- 
proximate 2-year  flood  in  March  1976  showed  the  flood  water 
to  be  clear  with  no  detritus. 
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With  the  preceding  information,  it  is  apparent  that 
Curve  1 of  Figure  11A  is  used  with  a D75  value  of  0.5  mm  to 
determine  the  non- scouring  velocity  for  all  of  the  samples 
with  no  plasticity  index.  From  the  chart,  a velocity  of 
39.6  cm/ s is  obtained.  This  value  must  be  corrected  for 
alignment,  bank  slope,  and  water  depth  (Figure  11B) . Align- 
ment and  bank  slope  are  assumed  to  be  constant  for  all  three 
discharges  under  consideration,  whereas  the  depth  of  flow 
increases  significantly  with  increasing  discharge.  For  the 
west  dike,  a bank  slope  of  4:1  gives  a slope  angle  (Z)  of 
~ 14°  with  a cotangent  of  : 4,  Therefore,  the  bank  slope 
correction  factor  is  not  applicable  to  either  the  dike  or 
the  floodplain.  The  alignment  factor  near  cross-section  8 
is  the  same  (assumed  so)  for  all  discharges.  The  curve 
radius  near  the  study  area  is  greater  than  1950.7  m and  the 
water  surface  width  is  = 121.9  m.  The  width  of  the  water 
surface  actually  does  change  with  discharge,  but  the  resul- 
tant change  in  the  alignment  factor  (curve  radius/water  sur- 
face width)  is  insignificant.  Depth  of  flow  for  the  dif- 
ferent floods  is  approximate  because  an  assumption,  that  the 
flood  channel  near  cross-section  8 is  rectangular,  is  made. 
This  assumption  is  valid  considering  the  small  contribution 
of  the  knoll  and  dike  side-slopes  to  the  total  area. 

The  water  depths  for  the  three  discharges  are: 

area  (m^)  width  (m)  depth  (m) 

2-year  flood  147.31  (147.54)  121.92  (121.92)  1.21  (1.21) 

10-year  flood  268.19  (267.39)  121.92  (121.92)  2.20  (2.19) 

100-year  flood  385.79  (383.88)  121.92  (121.92)  3.16  (3.15) 

NOTE:  DNR  setting  values  are  followed  by  Corps  of  En- 
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gineers  setting  values  in  parantheses. 

The  flood  channel  near  cross-section  8 is  slightly 
sinuous.  Hap  analysis  and  personal  observations  during  a 2- 
year  flood  lead  to  this  conclusion.  It  is  possible  that 
the  less  frequent  events  could  destroy  the  sinuosity  because 
of  the  greatly  increased  depths  and  volumes  of  water.  How- 
ever, for  the  purposes  of  this  research,  the  sinuosity  is 
assumed  to  remain  constant  for  all  events  under  considera- 
tion. 

With  the  data  heretofore  presented,  the  maximum  per- 
missible velocities  that  can  exist  before  scour  should  oc- 
cur are:  g 
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However,  these  velocities,  as  determined  by  the  SCS  method, 
are  not  sacrosanct.  The  agency  is  wary  of  conclusions  based 
upon  allowable  velocities  of  less  than  61  cm/s.  Personal 
communication  with  two  SCS  engineers  could  provide  no  ade- 
quate reasons  as  to  why  the  method  used  in  the  present  re- 
search should  not  be  valid  (Brown,  1976;  Hamper,  1976).  In 
order  to  get  another  view  of  the  problem,  the  critical  ve- 
locity method  based  upon  Russian  data  is  employed  (Chow, 
1959).  If  the  grain  size  value  is  used  with  the  Chow 
method  (Figure  12) , the  same  allowable  velocities  are  ob- 
tained as  with  the  SCS  method.  This  fact  is  not  surprising 
though,  since  the  SCS  method  is  partly  based  upon  the  same 
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Russian  data.  If  the  average  grain  size  (D^q  = 0.25  mm  or 
2 phi)  is  used,  different  velocities  are  obtained: 

TABLE  6:  Expected  Velocities  and  Critical  Velocities 

Expected  V Critical  V Critical  V 
(HEC-2)  (Chow)  (SCS) 

2-year  flood  18.3  cm/s  3 2 cm/ s 37.5  cm/s 

10-year  flood  27.4  cm/s  36  cm/ s 42.1  cm/s 

100-year  flood  36.6  cm/s  38.1  cm/s  44.5  cm/s 

Although  both  critical  velocity  methods  are  based  upon  em- 
pirical results,  this  researcher  is  confident  that  the 
critical  velocities  obtained,  the  SCS  values  and  the  lower 
limit  figures  from  the  method  described  by  Chow,  are  rea- 
sonable. The  critical  velocity  curve  of  Hjulstrom  (1935) 
cannot  be  used  for  comparison  because  it  is  based  upon  uni- 
form, monodispersed  material. 

There  appears  to  be  little  doubt  that  the  smaller  two 
floods  do  not  have  the  velocities  required  to  scour  the  west 
dike  and  the  sands  on  the  floodplain  in  the  vicinity  of 
cross-section  8,  no  matter  what  allowable  velocity  is  con- 
sidered. The  100-year  event  is  not  so  clear-cut,  however. 

The  difference  between  the  expected  velocity  and  the  allow- 
able velocity  obtained  with  the  Chow  data  is  too  small  to 
discount  the  possibility  of  scour. 

If  complete  faith  could  be  placed  in  the  SCS  allowable 
velocities,  then  it  would  be  safe  to  say  that  scour  would 
not  occur,  even  during  the  100-year  event.  But  the  doubts 
about  the  SCS  method  for  the  lower  velocities,  combined  with 
the  results  of  the  analysis  with  the  Russian  data,  require 
that  the  conclusion  to  this  part  of  the  hypothesis  test  re- 
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mains  uncertain.  Thus,  scour  could  possibly  occur  during 
the  100-year  flood. 

What  about  the  expected  velocity  values?  They  are  ob- 
tained from  the  HEC-2  analysis,  which  relies  upon  the  cur- 
rently accepted  discharge  values  for  each  recurrence  inter- 
val. If  the  climate  has  changed  toward  a wetter  and  cooler 
regime  since  about  1950,  it  is  possible  that  the  discharge 
values  used  in  the  HEC-2  analysis  are  too  low,  which  would 
mean  that  the  expected  velocities  are  probably  too  low  as 
well.  Before  these  matters  are  explored,  the  backwater  hy- 
pothesis must  be  tested  with  the  same  data  used  above  in  the 
scour  analysis. 


67 


CHAPTER  V 

THE  BACKWATER  HYPOTHESIS:  DATA  AND  ANALYSIS 

Backwater  is  an  important  consideration  in  any  large 
development  project  on  a river  floodplain.  Limitations  of- 
ten exist  as  to  how  much  backwater  can  be  created  by  a new 
construction.  A 15.24  cm  limit  exists  in  Wisconsin  for  the 
amount  of  increase  in  water  surface  height  permitted  near 
urban  areas.  Previous  maximum  values  for  the  backwater  ef- 
fect resulting  from  the  Columbia  Generating  Station  ranged 
from  7.9  cm  (Corps  of  Engineers)  to  15.54  cm  (Thomas  Lee  of 
DNR)  . 

The  backwater  question  at  Portage  is  significant  for  a 
number  of  reasons:  (1)  the  Columbia  Site  Impact  Study  Group 
wants  to  know  which  figures  are  correct,  and  why,  because 
their  research  objectives  include  the  requirement  for  a 
thorough  understanding  of  all  factors  to  be  considered  in 
the  siting  of  future  power  plants,  (2)  the  possibility  ex- 
ists that  enough  backwater  is  created  to  significantly  alter 
the  expected  stage-discharge  relationship  and,  thus,  affect 
the  expected  capability  for  flood  control  of  the  Portage 
levees,  and  (3)  the  need  exists  to  know  the  contribution  of 
this  project  to  the  water  surface  heights  at  Portage  so  that 
it  can  be  accounted  for  when  new  developments  are  to  take 
place  in  the  vicinity  of  the  study  area. 

The  results  of  the  backwater  investigation  (Tables  7A 
and  7B)  are  not  as  dramatic  as  expected.  Both  hydraulic 
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settings  (those  of  the  DNR  and  the  Corps  of  Engineers)  pro- 
vide maximum  water  surface  elevation  increases  far  less  than 
the  15.24  cm  calculated  by  Mr.  Lee.  In  fact,  the  maximum 
value,  6.4  cm  during  the  100-year  flood  at  cross-section  8.1 
in  the  Corps  of  Engineers  setting,  is  less  than  the  maximum 
value  espoused  by  the  Corps  during  the  1972  environmental 
impact  proceedings. 

The  reason  for  the  large  discrepancy  between  the  re- 
sults of  the  present  study  and  those  of  the  earlier  studies 
stems  from  the  interpretation  of  what  double  encroachment  is. 
Under  the  law,  double  encroachment  means  that  an  amount  of 
conveyance  equal  to  that  lost  on  one  side  of  a river  because 
of  an  encroachment  must  be  subtracted  from  the  other  side  of 
the  river.  In  other  words,  construction  which  is  identical 
to  that  proposed  must  be  assumed  to  exist  opposite  the  pro- 
posed plant  on  the  other  side  of  the  river.  The  backwater 
which  would  result  from  the  encroachment  of  both  plants  is 
that  which  must  be  considered  to  satisfy  the  terms  of  double 
encroachment  (Schmied,  1973).  Most  state  ordinances,  in- 
cluding that  of  Wisconsin,  consider  double  encroachment  only 
for  the  100-year  (regional)  flood. 

Results  of  this  study  depict  the  magnitude  of  the  ex- 
pected backwater  effect  for  three  floods  and  two  different 
hydraulic  settings.  For  the  DNR  setting,  double  encroach- 
ment is  accounted  for  in  the  study's  analysis  of  the  "Pres- 
ent" situation  (Table  7B) . The  significant  findings  are  two: 
the  magnitude  of  the  possible  backwater  effect  is  far  less 
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than  that  determined  by  Mr.  Lee  or  any  of  the  previous  re- 
searchers, and  the  greatest  effect  of  the  backwater  occurs 
near  the  power  plant,  as  the  effect  of  the  backwater  is 
greatly  diminished  near  Portage  itself. 

From  the  data  available,  it  seems  that  the  encroach- 
ment of  the  Columbia  plant  creates  a lesser  impact  than 
originally  thought.  The  ability  of  the  river  to  spread  out 
across  its  floodplain  during  flood  events  probably  accounts 
for  the  results  to  a great  extent.  The  small,  if  existent, 
backwater  values  for  cross-sections  16,  14,  and  13  near 
Portage  suggest  a rejection  of  the  hypothesis  that  the  en- 
croachment of  the  Columbia  plant  may  affect  the  Wisconsin- 
Fox  River  drainage  situation  during  floods  of  modest  recur- 
rence intervals. 

Finally,  it  seems  that  the  plant's  effect  during  the 
regional  flood  is  well  within  the  limits  established  by  the 
DNR  floodplain  zoning  ordinance.  The  reason  for  the  large 
discrepancy  between  the  results  of  this  study  and  others  in- 
volves the  interpretation  of  double  enroachment.  Only  the 
area  directly  across  the  river  from  the  Columbia  plant  was 
considered  in  the  double  encroachment  procedure  in  the  pres- 
ent study.  Earlier  research  done  by  Mr.  Lee  considered 
double  encroachment  from  the  power  plant  upstream  to  the 
city  of  Portage.  The  interpretation  of  the  ordinance  made 
in  the  present  study  is  believed  to  be  the  intent  of  the 
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The  6.4  cm  water  surface  elevation  increase  at  cross- 
section  8.1  (C  of  E setting)  is  a substantial  increase,  even 
if  it  falls  within  allowable  limits.  Whether  this  increase 
has  any  geomorphic  significance  is  questionable,  however. 

It  is  of  interest  to  consider  how  the  backwater  effect 
may,  in  fact,  be  of  substantial  importance  if  indiscriminate 
floodplain  development  occurs.  Multiple  developments  com- 
pound the  backwater  effect,  possibly  to  such  an  extent  that 
any  limits  are  greatly  exceeded  because  of  the  combined  ef- 
fects of  more  than  one  project.  The  possibility  of  such  an 
occurrence  in  the  study  area  appears  to  be  minimal  because 
the  floodplain  ordinance  is  shown  to  be  much  more  effective 
than  originally  thought. 


CHAPTER  VI 

THE  CLIMATE  CHANGE  HYPOTHESIS:  DATA  AND  ANALYSIS 

A change  to  a cold  and  wet  climate  regime  could  have  a 
significant  effect  on  flood  magnitude  and  frequency,  and 
thus  on  geomorphic  analyses  of  the  type  presented  in  this 
research.  Today's  climate  in  the  Upper  Wisconsin  River  Val- 
ley may  be  similar  to  that  of  the  cold,  wet  meridional  cli- 
matic regime  that  existed  prior  to  1895. 

1873-1910  DATA 

The  first  step  in  the  test  of  the  climate  change  hy- 
pothesis is  an  examination  of  rainfall  and  streamflow  data 
for  the  period  1873-1910  in  order  to  detect  the  change  in 
precipitation  regimes  thought  to  have  occurred  in  1895,  and 
also  to  determine  the  precipitation  characteristics  of  a 
meridional  regime. 

Monthly  rainfall  data  are  grouped  into  the  three  hy- 
drologic seasons.  Graphs  of  precipitation  versus  time  (Fig- 
ure 13)  depict  a noticeable  difference  in  precipitation 
values  between  the  first  and  second  halves  of  the  period  in 
each  of  the  three  seasons.  The  time  series  show  that  a 
shift  from  a wetter  to  a drier  regime  occurred  in  the  mid- 


1880s. 


Running  means  are  simple  ways  to  show  the  existence  of 


a change  in  climate  and  the  presence  of  a trend  in  climatic 
data.  The  5-year  (Figure  14)  and  10-year  (Figure  15)  run- 
ning means  depict  a significant  change  in  climate  from  a 


ure  13:  Precipitation  and  Fall  Low  Flows  1873-1910 


ear  Running  Means  1873-1910  Precipitation  and  Fall  Low  Flows 
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wetter  to  a drier  regime  in  each  hydrologic  season  (.001  > 

P for  Fall  and  Winter  and  . 01  > P for  Spring).  The  change 
appears  to  have  taken  place  around  1885  in  all  seasons. 
Graphs  depicting  cumulative  deviations  from  the  mean  (Fig- 
ure 16)  effectively  show  the  time  when  the  change  occurred. 
The  year  1885  appears  to  be  the  time  of  change  for  the 
Spring  and  Fall  hydrologic  seasons,  but  the  evidence  points 
to  1890  as  the  year  that  the  climate  changed  to  a drier  re- 
gime in  the  Winter  season. 

Streamflow  data  for  the  1873-1910  interval  are  limited 
to  low  flows  for  the  Fall  hydrologic  season.  Low  flow’s  are 
appropriate  because  they  should  accurately  reflect  the  de- 
gree of  wetness  or  dryness  of  the  prevailing  climate.  LowT 
flows  can  be  employed  because  the  Wisconsin  River  was  not 
regulated  to  any  extent  during  the  1873-1910  period.  Unfor- 
tunately, low  flow  data  for  only  the  Fall  season  are  avail- 
able. Too  many  low  flow  data  are  missing  fro.n  the  record 
to  use  the  low  flow  variable  for  the  Spring  and  Winter  sea- 
sons. However,  the  Fall  is  a season  of  great  significance 
because  most  of  the  high-magnitude  flood  events  on  the  Wis- 
consin River  have  occurred  during  that  season  (C  of  E,  1971) 
Low  flows  for  the  Fall  are  presented  in  the  same  fash- 
ion as  the  precipitation  data.  Graphs  of  low  flows  versus 
time  (Figure  13) , 5-year  and  10-year  running  means  (Figures 
14  and  15),  and  cumulative  deviations  from  the  mean  (Fig- 
ure 17)  are  employed  in  the  analysis.  Graphically,  Fall 
low  f 1 ows  correspond  quite  well  with  Fall  precipitation  data 
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(Figure  13) , although  a regression  analysis  produced  an 
of  only  +0.35,  .001  > P.  A graph  of  cumulative  deviations 
from  the  mean  (Figure  17)  clearly  shows  the  change  in  1885 
from  a wetter  to  a drier  regime. 

Although  the  data  do  not  reflect  a change  in  climate 
in  1895  as  proposed,  they  do  clearly  portray  the  existence 
of  a change  in  regime,  beginning  in  1885  for  the  Spring  and 
Fall  seasons  but  not  until  1890  for  the  Winter.  The  meridio- 
nal regime  is  much  wetter  than  the  zonal  regime  that  suc- 
ceeded it.  Persistence  of  variability,  a characteristic  Oi 
meridional  circulation  (Knox  et  al.,  1975),  is  most  notice- 
able during  the  pre-1885  time  of  meridional  circulation  in 
the  Spring  season.  The  Spring  meridional  regime  has  highly 
variable  precipitation  values,  all  above  the  mean.  However, 
the  Spring  variability  remains  persistent  throughout  the  en- 
tire period  of  record,  although  the  later  values  are  not  as 
great  as  the  pre-1885  precipitation  totals.  Persistence  of 
variability  is  not  as  marked  in  the  other  seasons,  as  there 
are  successions  of  years  with  persistent  wetness  or  dryness 
and  no  evident  persistence  of  variability.  Thus,  from  the 
rainfall  and  streamflow  data,  the  following  inferences  can 
be  drawn:  (1)  that  a shift  in  climate,  from  a wetter  to  a 
drier  regime,  did  occur;  (2)  that  the  change  in  climate  took 
place  around  1885,  rather  than  1895;  (3)  that  the  pre-1885 
regime,  the  meridional  circulation  regime,  is  characterized 


by  amounts  of  precipitation  in  all  three  hydrologic  seasons 
greater  than  the  totals  for  the  zonal  post-1885  period;  and 
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(4)  that  persistence  of  variability,  as  a characteristic  of 
a meridional  regime,  appears  to  be  less  significant  than  the 
increased  amounts  of  precipitation. 

1925-1975  DATA 

Lamb  (1966)  and  others  hypothesized  that  the  upper  air 
circulation  patterns  in  the  northern  hemisphere  shifted  from 
a zonal  back  to  a predominantly  meridional  regime  around 
1950.  Precipitation  and  streamflow  data  from  the  period 
1925-1975  are  analyzed  in  order  to  test  for  the  proposed 
change  in  climatic  regime.  The  hypothesized  meridional  re- 
gime should  resemble  the  meridional  regime  that  existed  in 
the  Upper  Wisconsin  River  Valley  prior  to  1885.  This  being 
so,  the  new  regime  should  have  greater  than  average  precipi- 
tation values  and  significant  persistence  of  variability, 
the  latter  characteristic  at  least  in  the  Spring. 

Precipitation  data  from  13  weather  stations  upstream  of 
the  Columbia  Generating  Station  (Table  1)  are  used  in  the 
analysis.  The  data  are  grouped  in  the  three  hydrologic  sea- 
sons. As  with  the  precipitation  data  from  the  1873-1910 
period,  these  data  are  also  refined  through  the  use  of  time 
series  (Figure  18A,  B,  and  C) , cumulative  deviations  from 
the  mean  (Figure  19),  and  5-year  and  10-year  running  means 
(Figures  20  and  21). 


An  examination  of  the  graphs  of  precipitation  versus 
time  for  each  hydrologic  season  reveals  no  conclusive  evi- 
dence of  a shift  in  precipitation  patterns  in  each  season. 


Figure  21;  10-year  Running  Means  1925-1975  Precipitation  and  Streamflow 
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However,  a trend  toward  a wetter  regime  since  about  1950 
does  exist  in  the  Fall  season.  Analysis  of  the  graphs 
which  depict  the  running  means  produces  the  same  results. 

In  fact,  it  appears  as  if  the  Winter  precipitation  has  de- 
creased since  1950.  Spring  precipitation  exhibits  year-to- 
year  variability,  but  this  variability  is  certainly  not  per- 
sistent. The  Fall  season  has  both  characteristics  of  a 
meridional  climate  regime:  greater  precipitation  and  per- 
sistent year-to-year  variability.  However,  the  persistence 
of  variability  carries  through  the  Fall  season  for  the  en- 
tire 1925-1975  period.  As  with  the  1873-1910  data,  the 
Spring  season  exhibits  the  best  evidence  of  persistent  vari- 
ability during  the  years  which  should  encompass  the  domi- 
nance of  the  meridional  regime,  approximately  1950-1975  in 
the  recent  period.  Winter  precipitation  is  again  anomalous, 
as  there  is  a trend  toward  low  precipitation  and  low  vari- 
ability since  1950.  The  same  information  can  be  gleaned 


from  graphs  of  cumulative  deviations  from  the  mean. 

Precipitation  data  for  1925-1975  offer  little  support 
for  the  hypothesis  that  the  climate  has  shifted  to  a meridio- 
nal regime,  characterized  by  higher  than  average  precipita- 
tion amounts  and  persistent,  high  year-to-year  variability. 
Only  the  Spring  season  exhibits  any  significant  persistence 
of  variability  that  can  be  dissociated  from  the  variability 
during  times  of  zonal  flow  domination.  All  three  seasons 
were  markedly  wetter  during  times  of  meridional  flow  in  the 


91 


terval  is  wetter  during  the  time  of  meridional  domination, 
1950-1975,  when  compared  to  the  zonal  circulation  period, 
1925-1949. 

Streamflow  data,  in  this  case  the  average  monthly  dis- 
charges from  12  gaging  stations  upstream  of  the  Columbia 
Generating  Station  (Table  2)  per  unit  area,  are  grouped  into 
the  three  hydrologic  seasons.  Graphs  of  flow  versus  time 
(Figure  18A,  B,  and  C) , cumulative  deviations  from  the  mean 
(Figure  22),  and  5-year  and  10-year  running  means  (Figures 
20  and  21)  are  employed  in  the  analysis.  Winter  flow  is 
sensibly  constant  throughout  the  entire  period  of  record, 
which  is  unfortunately  only  1941-1973  because  of  the  limited 
gaging  station  data  available.  However,  both  the  Fall  flows 
and  the  Spring  flows,  especially  the  latter,  exhibit  a trend 
toward  increasing  flow,  particularly  since  about  1960.  Be- 
cause of  river  regulation,  perhaps  most  emphasis  should  be 
placed  on  the  precipitation  data.  However,  the  manner  in 
which  the  streamflow  data  are  presented  is  the  most  accurate 
reflection  of  the  flow's  relationship  to  precipitation  and 
climate  regime  which  can  be  used  for  a river  that  contains 
dams  and  reservoirs.  The  streamflow  data  are  of  signifi- 
cance; their  only  real  shortcoming  is  the  length  of  their 
record. 

Analysis 


On  the  basis  of  the  precipitation  and  streamflow  data 
for  both  of  the  periods  analyzed,  1873-1910  and  1925-1975, 
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a number  of  statements  can  be  made:  (1)  that  the  change  from 
a zonal  to  a meridional  regime  around  1950  is  not  nearly  so 
well  reflected  in  the  data  as  was  the  opposite  change  in  re- 
gimes which  occurred  around  1885  (although  this  development 
may  be  influenced  by  some  bias  in  the  study's  period  of 
analysis  1873-1910  resulting  from  the  inclusion  of  the  de- 
cade 1875-1885,  a period  of  possible  exceptional  wetness  in 
the  Upper  Mississippi  Valley  region  (Knox  et  al. , 1975)); 

(2)  that  the  only  season  in  the  later  period  that  depicts 
both  of  the  characteristics  thought  to  accompany  a meridio- 
nal regime,  increased  precipitation  and  persistence  of  vari- 
ability, is  the  Fall,  and  even  then  the  results  are  less 
than  conclusive;  (3)  that  persistence  of  variability  carries 
through  both  periods  in  only  the  Spring  season,  but  there  is 
no  evidence  that  the  Spring  meridional  regime  is  wetter  than 
the  zonal  regime  in  the  1925-1975  period;  and  (4)  that  the 
Winter  season  is  anomalous  in  both  periods,  as  it  exhibits 
no  increased  precipitation  in  the  meridional  regimes  nor  any 
persistent  year-to-year  variability  that  can  be  dissociated 
from  the  variability  of  the  zonal  regimes. 

The  precipitation  and  streamflow  data  for  the  Upper 
Wisconsin  River  Valley  do  not  provide  conclusive  evidence 
for  a change  in  climate  around  1950.  Only  the  Fall  season 
shows  a trend  toward  above-average  precipitation  since  that 
time.  Persistence  of  variability  does  exist  in  the  Fall, 
but  it  occurs  throughout  the  entire  1925-1975  period.  Only 
the  Spring  season  has  variability  that  is  more  persistent  in 


94 


r 7 

the  meridional  regime  than  in  the  zonal  regime.  The  Winter 
season  has  no  characteristics  of  a meridional  regime.  Nei- 
ther greater  precipitation  values  nor  persistence  of  vari- 
ability exist  in  the  post-1950  Winter  data.  Any  variability 
that  is  evident  in  the  Winter  occurs  during  the  pre-1950  zo- 
nal regime. 

It  appears  that  the  Fall  season  is  the  only  time  which 
may  have  a significant  effect  on  hydrologic  data  because  of 
the  change  to  a wetter  and  colder  meridional  regime.  For 
the  study  area,  the  Fall  is  a season  of  great  significance 
because  the  dams  and  reservoirs  have  no  storage  effect  on 
flood  flows  during  that  time  (Sherrar,  1976a).  Since  the 
recent  data  exhibit  a trend  toward  increased  precipitation 
and  higher  discharges  in  the  Fall,  the  analysis  includes  a 
detailed  examination  of  the  rainfall  data  for  that  season. 
Precipitation  data  at  each  weather  station  (Figure  23)  for 
the  months  of  April  through  September  are  analyzed.  Before 
the  analysis  continues,  however,  a look  at  the  "normal"  pre- 
cipitation regime  for  the  Spring  and  Fall  seasons  is  appro- 
priate . 

Descriptions  of  the  "normal"  precipitation  regime  of 
Wisconsin  are  provided  by  Trewartha  (1961)  and  Barry  and 
Chorley  (1970).  Both  sources  depict  the  precipitation  rec- 
ord as  containing  a double  maximum  with  precipitation  peaks 
in  May- June  and  August-September . The  area  of  the  Upper 
Mississippi  River  has  a secondary  minimum  which  usually  oc- 
curs in  July-August  even  though  the  air  is  quite  warm  and 
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moist.  A ridge  of  high  pressure  over  the  Mississippi  Val- 
ley is  apparently  responsible  for  decreased  summer  thunder- 
storms and  the  resultant  secondary  precipitation  minimum. 
Normally,  cyclonic  activity  increases  again  in  early  Sep- 
tember as  the  polar  front  shifts  to  the  south  and  clashes 
with  warm,  moist  mT  (maritime  Tropical)  air  from  the  Gulf  of 
Mexico.  Frontal  wave  activity  increases,  making  the  first 
half  of  September  a rainy  period  in  Wisconsin.  After  about 
the  20th  of  the  month,  anticyclonic  conditions  return, 
bringing  warm,  dry  weather  conditions  commonly  referred  to 
as  an  Indian  summer. 

Trewartha  (1961)  concludes  that  the  summer  rainfall  in 
Wisconsin  is  primarily  the  result  of  large-scale  synoptic 
conditions  rather  than  the  consequence  of  random  thermal 
convection.  Barry  and  Chorley  (1970)  point  to  the  incursion 
of  polar  air  in  early  September  and  the  resultant  clash  with 
air  from  the  Gulf  as  the  cause  of  the  second  precipitation 
maximum  in  September,  although  Trewartha  (1961)  states  that 
in  northern  Wisconsin  the  same  maximum  occurs  in  August. 

The  secondary  minimum,  normally  in  July  or  August,  corre- 
sponds with  periods  of  great  anticyclonic  activity.  The 
Indian  summer  in  the  second  half  of  September  is  associated 
with  anticyclonic  conditions  resulting  from  strong  zonal 
flow. 

Bryson  and  Lahey  (1958)  hypothesized  that  the  sharp  in- 
crease in  precipitation  during  April,  which  culminates  in 
the  first  precipitation  maximum  in  May- June,  is  related  to 


a decline  in  the  zonal  index.  The  precipitation  maximum  is 
the  result  of  a change  in  the  wave  numbers  of  the  Rossby 
long  wave  patterns  with  an  accompanying  shift  in  the  posi- 
tions of  low  pressure  (trough)  locations.  Simultaneously, 
an  increase  in  frontal  activity  with  resultant  rainfall  oc- 
curs. They  attribute  the  mid-summer  precipitation  minimum 
to  the  northward  movement  of  the  Subtropical  High  Pressure 
cells  and  accompanying  prevalence  of  anticyclones.  In- 
creasing meridional  flow  in  early  September  causes  the  warm, 
moist  Gulf  air  to  rise  over  the  cold  cP  (continental  Polar) 
air  from  Canada.  The  results  are  the  high  precipitation  to- 
tals which  produce  the  second  precipitation  maximum  in 
either  early  September  or  late  August.  Increased  zonal  flow 
occurs  in  mid-September.  This  westerly  airflow  halts  the 
polar  air-Gulf  air  clashes  and  produces  the  warm,  dry  In- 
dian summer. 

Trewartha  (1961)  based  the  conclusions  of  his  study  on 
data  from  two  periods  of  time,  1906-1935  and  1920-1950. 

Both  periods  fall  within  the  pre-1950,  post-1895  zonal  cli- 
mate regime  proposed  by  Lamb  (1966)  and  Knox  and  others 
(1975).  Barry  and  Chorley  (1970)  do  not  refer  to  any  spe- 
cific data  base,  nor  do  they  cite  Trewartha  (1961)  or  Bryson 
and  Lahey  (1958)  in  their  bibliography.  Bryson  and  Lahey 
used  data  from  the  years  1899-1944  in  their  study.  These 
years  are  also  within  the  period  of  strong  zonal  airflow. 

The  descriptions  of  the  "normal"  summer  and  fall  precipita- 
tion characteristics  rely  on  alternating  periods  of  meridio- 
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nal  and  zonal  flow  to  account  for  the  double  maximum  and 
secondary  minimum  precipitation  characteristics.  When  one 
considers  that  the  data  bases  used  in  the  studies  were  all 
from  times  of  dominant  zonal  flow  (approximately  1885-1949), 
it  is  significant  to  consider  how  those  characteristics 
could  be  altered  if  the  data  included  precipitation  values 
fjLom  the  1950-1975  regime  of  meridional  circulation.  Since 
the  Fall  season  exhibits  a trend  toward  increased  precipita- 
tion and  since  during  the  Fall  the  dams  and  reservoirs  have 
no  storage  effect  on  the  flood  flow  of  the  Wisconsin  River, 
the  analysis  examines  the  precipitation  at  each  weather  sta- 
tion in  order  to  see  how  data  from  a meridional  regime  may 
affect  the  "normal"  precipitation  characteristics. 

Graphs  of  the  monthly  rainfall  from  each  of  the  13 
weather  stations  depict  how  the  zonal  data  from  1925-1949 
and  the  meridional  data  from  1950-1975  compare  (Figures  24A 
through  G) . The  zonal  flow  data  plot,  without  exception, 
in  the  manner  prescribed  by  Trewartha,  that  is,  with  a dou- 
ble maximum  and  a secondary  minimum.  All  stations  have  the 
first  peak  in  June  and  all  but  Rest  Lake  have  their  second 
peak  in  September.  All  stations  have  the  secondary  minimum 
in  either  July  or  August.  Thus,  the  precipitation  data 
match  the  "normal"  precipitation  characteristics  described 
above  at  all  stations  during  the  period  of  zonal  airflow. 
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However,  the  data  for  the  1950-1975  period  of  meridio- 
nal flow  possess  quite  different  characteristics.  The  most 
obvious  change  in  characteristics  is  probably  the  marked  re- 
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duction  in  the  magnitude  of  the  June  peak,  and,  in  some 
cases,  the  elimination  of  that  peak.  Another  significant 
change  is  the  shift  of  the  second  precipitation  maximum  from 
September  to  August  at  many  stations.  The  initial  peak  has 
also  shifted  from  June  to  May  at  some  of  the  stations.  With- 
out exception,  the  month  of  April  is  wetter  during  the  1950- 
1975  period.  A further  change  is  the  increased  precipita- 
tion during  either  July  or  August  or  both  at  all  of  the  sta- 
tions except  Coddington.  Finally,  precipitation  totals  de- 
creased in  September  at  all  but  two  of  the  stations. 

From  the  graphs  of  the  monthly  precipitation  at  each 
station,  it  is  quite  apparent  that  the  precipitation  charac- 
teristics no  longer  resemble  those  previously  considered  to 
be  "normal".  This  analysis  of  monthly  rainfall  demonstrates 
that,  indeed,  the  climate  in  the  study  area  has  changed 
since  1950.  A trend  toward  increased  precipitation  has  been 
shown  to  exist  in  the  Fall  hydrologic  season.  Monthly  data 
at  each  station  depict  that  the  timing  and  amount  of  pre- 
cipitation during  the  months  of  April  through  September 
have  changed  markedly.  The  characteristics  of  the  meridio- 
nal climate  of  the  post-1950  period  are  not  strikingly  simi- 
lar to  those  of  the  pre-1885  meridional  regime  (again,  this 
fact  may  result  partly  from  the  possibility  of  extreme  wet- 
ness during  the  1875-1885  period),  but  no  matter,  the  analy- 
sis has  demonstrated  without  qualification  that  the  charac- 
teristics of  the  Fall  and  most  of  the  Spring  hydrologic  sea- 
sons during  the  years  1950-1975  reflect  a different  climate 
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regime  from  the  characteristics  of  those  years  dominated  by 
strong  zonal  airflow. 
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CHAPTER  VII 
THE  REAPPRAISAL 

The  change  in  climate  since  1950  for  the  Upper  Wiscon- 
sin River  Valley  is  best  manifested  in  the  altered  timing 
of  "normal"  rainfall  during  the  months  of  April  through 
September.  The  Fall  hydrologic  season  has  recorded  higher 
precipitation  averages  since  that  time,  even  though  some 
months  in  that  season  have  not.  Some  months  during  the 
Spring  hydrologic  season  also  have  had  greater  precipitation 
averages  since  1950.  An  evaluation  of  the  impact  of  the 
demonstrated  change  in  climate  on  the  study  problems  is  ap- 
propriate . 

In  light  of  the  change  in  "normal"  precipitation  char- 
acteristics, the  discharge  values  obtained  from  the  USGS  and 
the  Corps  of  Engineers  may  not  be  representative.  The  in- 
tent of  this  chapter  is  to  develop  streamflow  data  for  the 
existing  climatic  regime.  Unfortunately,  no  simple  proce- 
dure for  the  proposed  task  exists.  The  rainfall  data, 
monthly  and  seasonal,  do  not  provide  a means  of  predicting 
discharge.  Normally,  storm  rainfalls  of  short  duration 
(hours)  or  steady  rainfalls  of  long  duration  (days)  are  em- 
ployed to  generate  peak  streamflow  values.  An  indirect  ap- 
proach to  the  problem  of  discharge  formulation  involves  a 
comparison  of  magnitudes  and  frequencies  for  two  time  peri- 
ods. A correction  factor  can  be  developed  to  account  for 
any  difference  between  the  two  distributions.  With  the  data 
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available  for  this  research,  only  the  indirect  approach  is 
feasible . 

The  data  for  the  analysis  are  the  maximum  gage  heights 
(stages)  of  the  Wisconsin  River  at  the  Portage  canal.  There, 
the  US  Weather  Bureau  maintains  a crest  gage  which  has  a pe- 
riod of  record  compatible  with  that  of  the  precipitation  data 
used  in  Chapter  VI.  Monthly  stages  for  the  1924-1975  period 
were  examined  in  order  to  produce  an  annual  flood  series. 

The  period  was  divided  in  half,  with  a zonal  interval  1924- 
1949  and  a meridional  interval  1950-1975.  Each  interval  was 
treated  as  a separate  population  with  its  own  distribution 
of  flood  frequencies.  The  formula  used  for  the  recurrence 
interval  was  r.i.  = N + 1/M.  Results  were  plotted  on  Gumbel 
probability  paper. 

Three  flood  frequency  determinations  were  made  both  for 
the  zonal  and  for  the  meridional  intervals:  one  for  the  Snow- 
melt period  (December  through  April) , one  for  the  Rainfall 
period  (May  through  November),  and  one  for  the  entire  year. 
The  Snowmelt  curves  (Figure  25)  show  larger  floods  for  most 
recurrence  intervals  since  1950.  Not  much  weight  can  be 
placed  on  these  results,  however.  River  regulation  and  ice 
jams  affect  the  Snowmelt  frequency  distributions.  Ice  can 
cause  backwater,  which  produces  misleading  stages  at  the 
crest  gage  (Dalrymple,  1960).  Reservoir  storage  is  a factor 
in  the  springtime,  as  the  reservoirs  are  carefully  monitored 
to  control  snowmelt  runoff  (Sherrar,  1976b).  The  yearly 
flood  frequency  curves  are  equally  suspect,  since  they  in- 
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elude  the  Snowmelt  data. 

The  flood  frequency  determination  for  the  Rainfall  in- 
terval is  not  affected  by  either  of  the  aforementioned  prob- 
lems; therefore,  it  should  provide  a reasonable  reflection 
of  the  climate  during  1924-1975.  With  the  a priori  knowledge 
gained  from  the  climate  change  analysis,  one  would  expect  the 
flood  frequency  distributions  for  the  zonal  and  meridional 
periods  to  differ  from  one  another.  A graphical  comparison 
of  the  two  distributions  (Figure  25)  shows  virtually  no  dif- 
ference between  the  two  curves.  To  check  the  situation  fur- 
ther, all  events  with  stages  greater  than  flood  stage  (5.2  m) 
for  the  interval  1909-1923  were  incorporated  into  the  zonal 
regime  data  by  a method  that  considers  historical  floods 
(Benson,  1962).  No  data  are  available  for  1900-1908.  Data 
from  the  years  before  1900  are  not  valid  in  this  analysis, 
because  the  Portage  levees  in  their  present  form  did  not  ex- 
ist until  1900  (EIS,  1974).  Again,  a graphical  comparison 
between  the  zonal  and  meridional  distributions  shows  virtu- 
ally no  difference  (Figure  26) . 

The  Portage  stage  data  indicate  that  the  meridional  re- 
gime has  had  no  noticeable  effect  on  the  flow  of  the  Wiscon- 
sin River  in  the  study  area.  Two  possibilities  exist:  either 
the  new  climate  has  had  no  effect,  or  it  has  had  an  effect 
which  cannot  be  determined  with  the  method  and  data  used.  In 
either  case,  representative  discharges  for  the  post-1950  cli- 
matic regime,  that  differ  quantitatively  from  those  of  the 
pre-1950  regime,  cannot  be  determined  in  the  course  of  the 


Figure  26:  Flood  Frequency 
Distributions  Extended  for  Portage 
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present  research.  This  is  unfortunate,  because  the  graphs 
of  the  monthly  precipitation  at  each  weather  station  (Fig- 
ures 24A-G)  indicate  antecedent  conditions  for  the  meridional 
regime  that  are  markedly  different  from  those  of  the  zonal 
regime.  For  geomorphic  problems,  such  as  scour,  antecedent 
conditions  are  primary  considerations . Because  of  the  fail- 
ure of  the  quantitative  effort  at  discharge  prediction,  only 
subjective  thoughts  as  to  the  effects  of  the  climate  change 
can  be  presented. 

Some  of  the  changes  in  precipitation  characteristics 
since  1950  are  applicable  to  all  stations:  Junes  and  Septem- 
bers are  drier  than  before,  while  Aprils  are  wetter  at  all 
but  two  stations.  In  the  northern  portion  of  the  Upper  Wis- 
consin Valley,  Augusts  are  wetter  at  all  stations  (Figure  23, 
stations  1-6)  and  Julys  are  wetter  at  three  stations,  drier 
at  two,  and  unchanged  at  one.  The  situation  in  the  southern 
portion  (Figure  23,  stations  7-13)  differs  somewhat,  as  Au- 
gusts are  wetter  at  only  three  stations,  but  Julys  are  wetter 
at  all  but  one.  Precipitation  averages  during  May  have  re- 
mained constant  throughout  the  entire  1925-1975  period  at  all 
stations . 

Changes  in  the  timing  of  "normal"  rainfall  should  have 
had  effects  on  many  geomorphic  processes  in  the  study  area. 
The  indicated  changes  could  have  affected  thresholds  of  sta- 
bility, and  thus  have  had  a major  impact  on  the  results  of 
any  analysis  of  the  possibility  of  scour.  Soil  moisture  con- 
ditions could  have  been  so  altered  that  discharges  that  were 
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once  geomorphically  ineffective  would  now  have  the  ability 
to  do  much  geomorphic  work.  The  changed  precipitation  re- 
gime, with  its  different  antecedent  co  ^icions,  might  cause 
more  frequent  high-magnitude  floods  during  the  Rainfall  in- 
terval than  the  present  flood  frequency  determinations  would 
predict.  A comparison  between  the  zonal  and  meridional  pe- 
riods of  the  number  of  days  when  the  Wisconsin  River  was 
above  flood  stage  (flood-days)  at  Portage  is  worthy  of  con- 
sideration in  this  regard. 

An  examination  of  the  number  of  flood-days,  obtained 
from  the  US  Weather  Bureau's  Daily  River  Stages , provides  in- 
conclusive results.  For  the  zonal  interval  1909-1949,  there 
were  40  flood-days,  exceeding  by  9 the  31  flood-days  of  the 
1950-1976  meridional  period.  Because  the  length  of  record 
for  the  two  periods  differs,  it  may  not  be  apparent  that,  on 
the  average,  most  flood-days  per  Rainfall  interval  have  oc- 
curred during  the  meridional  climate  regime.  However,  the 
difference  in  the  averages  between  the  two  periods  is  small-- 
1.15  flood-days  per  Rainfall  interval  in  the  meridional  re- 
gime compared  to  0.98  flood-days  per  Rainfall  interval  in  the 
zonal  regime.  An  analysis  of  the  temporal  distribution  of 
flood-days  is  of  greater  value  (Figure  27) . A readily  ap- 
parent observation  is  the  striking  increase  in  flood-days 
during  May,  from  none  in  the  zonal  period  to  20  in  the  me- 
ridional interval.  On  the  other  hand,  the  number  of  flood- 
days  in  June,  September,  and  October  decreased  in  the  meridio- 
nal period.  Importantly,  the  changes  in  the  timing  of  flood- 
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Figure  27:  Temporal  Distribution  of  Flood  Days  at  Portage 
(for  months  of  April  through  November) 
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days  at  Portage  corresponds  quite  well  with  the  various 
changes  in  the  timing  of  precipitation  maxima  and  minima 
in  the  Upper  Wisconsin  River  Valley  (See  Figures  24A-G  and 
27).  It  is  surprising  to  note  the  absence  of  any  flood-days 
in  August,  especially  considering  the  precipitation  increases 
throughout  the  Wisconsin  River  drainage  basin  above  Portage 
during  that  month. 

If  the  April  flood-days  could  be  included  in  the  analy- 
sis, the  difference  in  the  number  of  flood-days  between  the 
two  climatic  regimes  would  be  noteworthy.  In  fact,  there 
would  be  an  average  of  2.51  flood-days  per  Rainfall  interval 
in  the  meridional  regime  and  only  1.65  events  per  Rainfall 
interval  in  the  zonal  period.  However,  April  is  a month 
which  may  have  floods  resulting  from  either  snowmelt  runoff, 
excess  rainfall,  or  a combination  of  the  two  preceding  causes. 
Because  of  the  uncertainty  about  the  origin  of  the  April 
flood  stages,  the  data  from  that  month  cannot  be  included  in 
the  analysis. 

The  small  difference  in  the  number  of  floods  per  Rain- 
fall interval  between  the  two  time  periods  provides  no  sup- 
port for  the  hypothesized  changes  in  hydrologic  response  ex- 
pected with  a shift  to  a meridional  climate  regime.  The 
changes  in  the  temporal  distribution  of  flood-days  may  pro- 
vide such  support.  It  appears,  in  the  Upper  Wisconsin  River 
Valley  at  least,  that  the  climatic  change  of  1950  is  mani- 
fested not  only  in  the  changed  timing  of  "normal"  rainfall, 
but  also  in  the  altered  timing  of  flood-days.  Such  a devel- 


opment  emphasizes  the  necessary  consideration  of  antecedent 
conditions  in  the  study  of  geomorphic  processes. 

The  analysis  in  Chapter  IV  of  the  scour  potential  dem- 
onstrated that  a possibility  exists  for  the  occurrence  of 
scour  during  the  100-year  flood.  The  discharge  value  of 
2690.4  m^/s  used  for  the  100-year  flood  may  be  too  small. 

The  Corps  of  Engineers  had  recognized  a value  of  2973.6  m^/s 
before  agreement  with  the  USGS  on  the  smaller  value  (EIS, 
1974).  A discharge,  not  much  larger  than  the  one  used  in 
the  present  study  for  the  100-year  flood,  could  produce  the 
threshold  velocity  necessary  to  scour  the  sandy,  unvegetated 
areas  of  the  floodplain  and  west  dike.  Computer  runs  for 
two  large,  hypothetical  discharges  provide  some  idea  of  ex- 
treme velocity-discharge  relationships  at  cross-section  8. 
Velocities  of  48.8  cm/s  and  54.9  cm/ s resulted  from  discharge 
values  of  4493.0  m~V s and  5380.8  m^/s  respectively.  Some- 
where between  the  lowest  of  these  discharge  values  and  the 
value  of  the  100-year  flood  used  in  the  present  research  lies 
that  discharge  which  produces  just  the  velocity  required  to 
initiate  particle  motion  in  the  study  area.  However,  at 
this  point  everything  is  conjecture.  A rigorous  analysis  of 
existing  data  leads  to  no  conclusive  evidence  for  the  occur- 
rence of  scour  or  excessive  backwater,  even  when  the  climate 
change  since  1950  is  taken  into  account. 
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CHAPTER  VIII 
CONCLUSIONS 

Scour  along  the  west  dike  and  on  the  floodplain  sur- 
face near  cross-section  8 may  occur  during  the  100-year 
flood.  There  appears  to  be  no  possibility  of  scour  during 
the  2-year  and  10-year  events.  The  backwater  effect,  as  a 
result  of  the  encroachment  of  the  Columbia  power  plant, 
does  not  exceed  the  limits  established  by  the  DNR  floodplain 
zoning  ordinance.  The  greatest  increase  in  water  surface 
elevation  is  6.4  cm.  This  value  occurs  at  a cross-section 
below  the  city  of  Portage,  as  the  backwater  effect  dimin- 
ishes upstream  from  the  encroachment.  The  backwater  effect 
has  no  significant  impact  on  the  normal  stage-discharge  re- 
lationships at  Portage  or  on  the  Wisconsin  River-Fox  River 
drainage  situation  during  flood  flows. 

Bogardi  (1968)  recognizes  four  possible  mechanisms  for 
the  failure  of  a dike:  wave  action,  seepage,  saturation, 
and  scour.  All  four  apply  to  the  study  site.  Waves  lap 
the  dike  on  windy  days.  Water  seeps  from  the  cooling  lake 
through  the  dike.  Floods  on  the  Wisconsin  River  are  of  long 
duration  (days).  However,  the  present  research  examines 
only  the  last  mechanism- -scour  from  the  effect  of  water 
flowing  against  the  dike  and  on  the  floodplain  surfaces. 

The  analysis  of  scour  potential  employs  the  assumption  that 
the  floodplain  and  west  dike  are  unvegetated.  In  reality, 
they  are  vegetated.  New  methods  to  determine  accurate 
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roughness  coefficients  for  vegetated  channels  are  being  in- 
vestigated by  Shen  (1973).  The  results  of  his  work  may  per- 
mit the  elimination  of  assumptions,  such  as  the  one  aforemen- 
tioned, and  permit  researchers  to  investigate  the  possibili- 
ty of  scour  with  the  vegetation  taken  into  account. 

A possibility  of  vegetational  change  exists  in  the 
sedge  meadow,  because  of  thermal  and  other  effects  resulting 
from  the  power  plant  (Willard,  1977).  The  occurrence  of  such 
change  could  affect  the  flow  hydraulics  of  floodwaters  near 
cross-section  8.  Floating  mats  of  peat  during  high  flows,  a 
troublesome  possibility,  could  also  affect  expected  velocity 
distributions  in  the  study  area.  Erosion,  as  a result  of 
overland  flow,  is  possible  on  the  west  dike.  The  dike  soils 
have  high  erodibility  indices  (Wischmeier,  1973),  but  the 
present  study  does  not  consider  the  effects  of  overland  flow. 
Future  research  is  needed  to  monitor  each  of  the  aforemen- 
tioned possibilities. 

The  backwater  question  is  relatively  straightforward. 

An  area  of  confusion  is  the  definition  of  double  encroach- 


ment. Double  encroachment  requires  that  the  amount  of  con- 
veyance lost  because  of  the  encroachment  be  subtracted  from 
the  total  amount  of  water  conveyed  on  the  floodplain  oppo- 
site the  encroachment.  In  other  words,  the  effect  of  one 
floodplain  project  is  determined  by  calculating  the  effect 
of  a similar  project  on  the  opposite  side  of  the  river.  The 
combined  effect  must  be  considered  (Schmied,  1973).  Pre- 
vious studies  (EIS,  1974)  considered  double  encroachment 
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along  the  entire  reach  of  the  study  area,  from  the  power 
plant  to  the  city  of  Portage.  The  present  study  considers 
only  the  area  directly  across  from  the  power  plant.  The  in- 
terpretation of  the  backwater  question  presented  herein  fol- 
lows the  intent  of  the  law. 

The  analysis  of  climate  change  produces  interesting  re- 
sults, but  the  results  do  not  lend  themselves  to  easy  inter- 
pretation. The  post-1950  meridional  regime  is  not  as  dra- 
matic in  character  as  the  pre-1885  regime.  Only  the  Fall 
hydrologic  season  exhibits  a trend  toward  increased  precipi- 
tation since  1950.  The  characteristic  of  high  year-to-year 
variability  is  difficult  to  dissociate  from  the  pre-1950 
data  in  all  three  hydrologic  seasons.  Changes  in  the  timing 
of  "normal"  precipitation  during  the  months  of  April  through 
September  are  amply  demonstrated.  The  different  antecedent 
conditions  should  affect  the  stream  discharge  during  those 
months.  Whether  the  change  in  antecedent  conditions  is  suf- 
ficient to  cause  a biogeomorphic  response  (Knox,  1972),  even 
on  the  smallest  scale,  is  questionable.  Over  an  extended 
period  of  time,  changes  in  antecedent  conditions  are  proba- 
bly more  likely  to  have  a geomorphic  effect  in  arid  areas 
as  opposed  to  humid  areas,  such  as  the  Upper  Wisconsin  River 
Valley.  Even  so,  such  changes  should  affect  geomorphic 
thresholds  in  the  study  area  and,  in  turn,  the  complex  re- 
sponses of  geomorphic  systems  (Schumm,  1973) . Further  re- 
search is  required  to  determine  the  extent  of  such  develop- 
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The  present  study  is  an  attempt  to  determine  the  abili- 
ty of  a stream  during  flood  flow  to  do  geomorphic  work.  The 
study  develops  a method,  albeit  in  need  of  refinement,  of 
the  type  called  for  by  Dury  (1973).  The  method  permits 
evaluation  of  the  potential  of  a stream  to  scour  its  flood- 
plain  and  floodplain  structures  during  floods  of  various 
magnitudes  and  frequencies.  Such  a method  should  be  useful 
in  fluvial  geomorphology,  let  alone  in  a school  of  neocatas- 
trophism.  The  analysis  of  climatic  data  should  provide  a 
basis  for  further  research  into  the  effects  of  climate 
change  in  the  Upper  Wisconsin  River  Valley.  Most  of  the 
goals  of  this  thesis  were  not  attained  without  qualifica- 
tion, which  further  demonstrates  the  complexities  involved 
in  the  treatment  of  geomorphic  processes.  Hopefully,  future 
developments  in  geomorphology  will  allow  similar  studies  to 
achieve  their  ends  without  qualification. 
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APPENDIX  1 

Particle  Size  Analysis  (Sieving) 
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NOTE:  All  Weights  in  Grams 
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APPENDIX  3:  continued 
Liquid  Limit  Determination 
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NOTE:  Results  with  Expedient  Method. 
All  Weights  in  Grams. 
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APPENDIX  4:  Precipitation  Data  1873-1910 


d r'lnotNr^oajncohvr^orv-^vooo'ONinNnMucMON'jm 
i ctf  <j-inr^vovo<T'Ocrir^voiDvDcocsJOr^vD<j-ni— icnictnooooc^ooooCT' 

O 0)  

i — i JS  st  -j-  -4-  st  u">  -tf-  sj-  vj-  <t  -4-  -<f  cn  cn  cn  co  co  cn  cm  cm  cm  cm  cn  cm  cm 

PS 

w • 

H d d I HOcOONfsfO'^OoOHr-KNOinNN 

f3  >*,  «J  I I CS|0\4sf<t\DH<JicOO\HOoOI^<tOHoONMnHOCT'^OCr\ 

M I <U  II 

5ir>  S • l <J-cl'4''4‘-<t'd-in'^r<l-'4‘vnLri^<l-'d-'<}'<t<-orOfOc^rorocNjCNlcS 


HrHX<X)Hr-IOOnOstvDvOHO>^OfOin^O\0(0^iDN-jN<f 
O inHmaiOOO\vOOMTv04'0  040\OOONr-IOONMO'vOCM^^O\D 

> 

<;  (’Kt't'tni-i'O'^'j'iri'J'iONt'ininfOincO'j-n^fnNcoNfNiNN 


< 

$ 


cn  ■J' m VO  rv  CO  o\  o I- ICMrO'J-in  cOr^KK^OHMCOsf-mcOI^COO'O 

r^i^iv.Nrsr^r^cooooO(Xioo  oooooooocoo'Cmtv<^omj'0'0\o\onO 
00  ON 


APPENDIX  4: 


1 


T3 

O 

0 

a 

m 

u 

C 

o 

o 


> 

o 


VI 

o 

S3  i cO 
MOO 
pShS 
fv, 

CO 


l O 

o o 
M X 


Pi 

w • 
HUC 
S3  >s  « 

M I O 


Vi  C NcorvsfOcnNvj  i i 

nj  iO  N<t  h <t  n <f  i i 

I O II 

u~>  s OHO'O'O'OOMN  I I 

H H 

mOO^VOr^OCNJr-lrH 

O coinN^cor^rocoooo 

^ 

(^vOCNOOOr^CT'VOvON 

1 — I r— I 

• • ooono>CT\cncjooioo 

S > NH(OiO4imOC0>J-H 

3 o 

U P>  HrlrlN^lOn^n  I 


ON  I I I I I I I I I 

'O'  I I I l I l I l l 

• I I I I I I I I I 

■ — *11111111 


V-I  p CM  lO  00  i — llOlOCsIOO  I I 

>s  10  CO  M M lO  CSI  CO  cr>  00  I I 

I O II 

m S M M 00  00  00  00  M lO  I I 

O vo  n o ouo  r-'  oo  m oo  m 

> 

C lO  M 00  00  O'  OO  LO  00  ID  'O' 

• • i£>  <T>  i— I C7>  CO  O' 'O' LO  CM  O 

6 > CIvO'OO'OiOHHHN 

3 O 

OQ  lO  lO  vt  CM  CM  CM  M i— I I 


m i i i i i i i i 

o>  I I I I I I I I 

• i i i i i i 

CM  i i i i i i 


i 


i i 
i i 


i 


CM  CM  CM  rn  <N  CM  O CO  I | 


CM  CM  N csl  Oi  CO  4 OO  00  1TI 


129 


u « 

• 

Vi 

Mnc^ooooooinroH 

rH 

SsC 

vO 

■ i i i i 

1 l l l 

i nj 

• i i i i 

l l 1 1 

o o 

• 

• i i i i 

1 1 l 1 

FALL 

1 

M 

oo 

• i i i i 

1 l l l 

APPENDIX  4:  Precipitation  Data  1925-1975 


s 

o 


►J 


3 (U 
U Q 


Po  C 
i nJ 
O 0) 

a 


e c 

l <D 
to  2 


O 

< 


§ > 
3 QJ 

C_>  Q 


o >,c 

2 I 
HO  0) 
Pi  r-t  S 
(X, 

C/3 


u a 

>s  ctf 
I (U 
to  £ 


O 


§ > 
3 a) 

C_>  Q 


>-l 

>s« 
I H) 
O <U 
r-t  a 

pr{ 

w • 

HUC 
Z !M« 
H I <U 

^ to  a 


130 

rH  0'3’l^'CNlDv£)tr)OcOCNiOiOrHCMvOCT\0»OiO<Ti>d  lOOOrOCSICy>COCO 
St  HN'JNOvOcOSO<tHOOOOOONH(NHHN  IOOOOOHO 

I (ON  StOrOOHH  rH  CO  pH  st  rs  vO  CO  CSI  iH  CO  CO  to  CO  to  OO 

I I I I I I I I I 

vO^vOOOtDcOOO'd'MON'J  vO\OcOCOiONCSI^I^iOOlOOO'd'h'0'0' 
NON'J-vf(OCO'^OOCOiOOOOCOr^O<SHO\  lOCO  OOCTNCTNOOCNOOOO'irH 

oooooooooooooooocr>CTvCT'oocooocooooorsrsr'»rsr'-rsr-^ooooovooav 

I I r- liOvOOOtONOi/KOOCMLDO't  OCM>JOO'ON>J  lOiOOcOOt 

I 

II  

i i OiO'cocotor^coooooaoCT'cocT'OOM^a'OONvor^Nr^r'r^coco 

i — i 

LOrsHOlHONiONO'vOO'vOtMtOOOcO  OvCSJ  O N N CO  N OO  ffi  N is 
COCSJstStO.HCO'3'l lOfsCTi  LO  sf  H Ost  LOvt  CO 

00(NOOtOrsLOO\LOiDOO|soONvO|s(NHNO|sS|siriaiLOO|siO 
i — I i I i — I i — I t — I i — I i — I H 

LOIsinCN-^  vOvOMOOCOHvOcO'J'JOON-^OiOLDffiCO'J-iO'd'OOlO 
r-HOO'IsrsMO'ONOO'O’LOStOOOOOHCO'J  NCMOOtOHst  OO 

HH(0(MH(OiOLO|srsONcOrssJ-IOCM  r- I Csl  CN|  1 i— I r-l 

llllllllllll — ll — I 

I I I I I I I I I III 

H\OHCOvO(^'4’r^COcOCO<tsfNHLOOH  COOOO  CM  ON  0<N  OOON 

aioicotOHr-isfioaiN'O-voaiH  ioiooco  oonh  aiootois\ooO'0’'i- 
\o\0'0'Oisis|sNisooooooooo'oooococoooisoorsrsisrsrsisrsN 

I I CTvstOVCOCTVCOOLncOstrs.rHavin<S|CTVvO<y«iH£SJstvOststOstCT\ 

I I VO  lO  CS  H COH  IsoOiOHvfcOrHLOLO 

II  

i i rs.r^rs.Lors.vovovOLorsr^oocxiCT'CT'C^ooc^coO'rsrs.r^rsoooooo 

st  ctics4<Nooooo''tfrHa'CNvOstcocnovOLOoorostr-icsioo  CM^ai-trs 
OrH'd-OCNCSIOOrsCSl(T>vOvO'X>cOtOoOO(T>  cost  HCOIOHOHHOO 

|sr^|sLOCOOOtOiOMO|siOiOHOOOIsO\  CTxOO  Ov  00  rs  lO  rs  |s  O N Oi 

i — 1 rH  rH 

i vo  vo  to  oo  rs  cm  lo  oavcorsLOooavLOvoovOcooostcovor^avst 
I CO  CO  i — I CO  St  Ov  to  «tOiOCOrsMNHcO  (OIS  HCOOOCOOuOOSOV 

I 

I rH  CS  OrHH  rHiHCOCNlCOCNCsIrHCNIststStstLOvOLOLO 

II  till  I 

I rHiHOOvOvOCOvOCOcOrHiHstCOvOcOOOoO  COvO  COOOOCMONOOMs 
I lOintOtOrsrsCHOvOMOCOSvOinisvOMsMnHHHOOvNCO 

I CO  CO  CO  CO  CO  CO  st  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  COCO  CO  CO  CO  CO  CO  CO  CS|  C\|  CS| 

I I I OOLOCOOrH-HCOvOStcOrHrHCOOOOCOCTvOOvOiHvOrHCOrHlO 
I I IrHCT'COCOCSICTvHrsrscOvOOVvOiOrHcOlscOfsOrHvOvOI'^iOO 

I II 

I I I CO  CN  CO  CO  CO  CO  St  CO  CO  CO  CO  CO  CO  CO  COCO  CO  CO  st  St  CO  CO  CO  CO  CO 

I iTiOrtfsrlvOlsvOHOiOO'OOcOrlNN  COCN  ■tlsOiOstO'HNvO 
i ostavoocoavr-^criOoocrirHrsCTvstcMvjDstr^ooostavoooorHio 

I * 

I COcOcOCNC^iHLOcOCSlLOCslstcOstcSstCsl  COCH  COvOcOCSlststcOcOCO 

invOlsoOO'OHMCO'tinvD|scOO'OH(S|CO«tirivONODONOHCv|(0 

CMCMCSiCNCNICOCOCOCOCOCOCOCOCOCOststststStstStStStStLOLrimLO 

CTv 


APPENDIX  4: 


131 


§ > 
3 <u 
o Q 


Vi 

>>C 

i cd 

o <u 
rJrH  g 


Vl  fJ 
So  cd 
i <u 
m g 


t3 

a) 

3 

a 

•H 

4J 

c 

o 

o 


o 

< 


B > 
3 <u 
u a 


O So  fS 
£3  i cd 

M O <U 
pirH  g 
Ph 
CO 


Vi  c 

>n  cd 
i a) 
m g 


O 

3 


6 > 
3 a) 
o Q 


Vi 

>n  C 
i cd 
o <u 

rH  g 

Prf 

Cd  • 

H Vi  fS 
!3  >v  cd 
M i CU 
^ m g 


o 

> 

c 


3 

w 

>1 


vD<j-'4-^0r-|cr>cM00i—ir^Lri<j-r^0NCN<}-o<j-c30r^c0'X> 

vor^cni-ih-r^oocoa\ff>a\OcocMcMs}-c£)OiAcoN<j 

m m r-^  on  oo  cm  cm  rH  NNHHNHNHr-incMH 

i i i i i i i i i i i i i i i i i i 

r-~  <f  CM  |-~  to  O On  on  <t  CM  u~>  U0  CM  1 I I I I 1 I I I 

cm  O cm  co  ^d- m oo  oo  cm  co  cm  on  i i i i i i i i i 

1 I 1 I I 1 I I I 

ON  On  ON  ON  ON  ON  00  CO  CO  ON  On  ON  CO  I 1 I I I I I I 1 

COOO'JmiOCOHHMCMHCONvOvOM'COCOOS  I I 
O ON  VO  CO  CO  ON  <*  ON  I",  r-l  CO  Ml-  r-l  CO  NO  CO  ON  m VO  r-l  I I 

II 

CO  1^- CO  ON  ON  ON  O ON  CO  ON  CO  CO  ON  CO  CO  CO  ON  ON  on  ON  I I 
i— I 

000NN0«JCMC0M,HCMrH0NM»1inc0MOCN|0'incM0N 

HNOONHNcDNCTNinr^t^cD4cOCOiClNOcOCOONOst 

HcoONONM’oooNONNOcoocONa'rNONONcooor^rN 
rH  i — I i — I I — I 

HNOCOCOiOUNCO'JOONCOMO'rNONCOONcOONinONCO 
ON  in  ON  rH  U1  t— I CO  nJ  CO  Cl  1C1  O CO  ON  NO  r- IvOCOUOHCO 

CO  CO  CO  r— I rH  iH  CM  CO  CM  CM  rH  rH  CM  rH  rH 

I I I I I I I I I I II 

NO  VO  VO  CO  ON  O CO  CM  CM  ON  CO  in  CO  I I I I I I I I I 

cm  on  o on  rH  oo  o co  m o rH  on  i i i i i ■ i i i 

I I I I I I I I I 

NiOrwONNOONNlNCOCOIN-  I I I I I I I I I 

in  O VO  0O  CO  rH  OO  <t  Ml- rH  <t  ON  00  rH  CO  00  O O ON  vt  | | 

rH  o cm  m on  oo  on  cm  co  cm  o co  co  co  oo  oo  in  in  i i 

oooorwovovovorMrviMrMNcooooocoiMiMrNrx  i i 

inMrN.inHOcMooNjHHcoincMvoooHvtiMcoMoo 

NncncMMcoHONHNiMinHONinnMoooNvonco 

ON|M[^vOinvOONvONvOvOCOiOCOHOOIMininHSvO 

I — I r — I 

CNHONvoinM-HinHiorNooiMcocoinNCDtOvl  co 

NfOMncovooovOdinM'NcOHsfO'jininoiM  -o 

lO  CO  CM  CM  H r-H  CM  CM  rH  CM  rH  CM  O 

I I I I I I I I I I 

rvi — d- cm  cm  cm  oo  co  in  co  rH  rH  i i i i i i • i i 

vOlONOHOCMOOtOvOaOIN  I I I I I I I I 1 

I I I I I I I I I 

CM  CM  CM  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  I I I I I I I I I 

CM  CM  -d-  VO  -^-  Ml-  CM  ON  ON  -O  co  CO  OO  VO  CO  CO  rH  VO  00  H I I 

MNin'dTNOoosvorvcvicMNvONdininoN'OO  i i 

II 

CM  CM  CM  CM  CM  CM  CO  CM  CM  CM  CO  CO  CO  CO  CO  CO  CO  CO  CO 'd-  I I 

incooocooNONoo.j-NinoO'd'O'dooNONoovovoiV'd' 

ONOOOHvOinHinHvtOHCOHOOONCNICOOOvOOO 

CMrHCMCOCMCMCOCM'd-rHCM'd'tH'd'rHcOCMin'd-'d'CMcO 


MinvONoooNOHcMrO'dinvorNoooNOHCMcoM'in 
minminminvovovovovovovovovovon-r^i — r-  r-~ 

ON  ON 


132 


APPENDIX  5:  Streamflow  Data  1873-1910  (in  m) 
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APPENDIX  6:  HEC-2  Data 

All  of  the  HEC-2  data,  to  include  both  computer  input 
and  output,  can  be  found  at  the  office  of  DNR's  Floodplain- 
Shoreland  Management  Section,  3rd  floor  of  the  Pyare  Square 
Building  on  University  Avenue  in  Madison,  Wisconsin. 
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